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Abstract 

The study of photoproduction reactions has gained a unique opportunity with 
the HERA data. The high center of mass energy allows to carry out quantitative 
tests of QCD, to explore the substructure of both the photon and the proton and to 
shed new light on soft and diffractive processes. In this report we review the HERA 
results, comparing them with the reach of fixed target photoproduction experiments 
and complementary measurements from e"*'e“ and hadron colliders. 


1 Introduction 


High energy photoproduction processes occur with high cross sections at the HERA ep 
collider (a/s = 300 GeV), due to the large flux of quasi real photons from the electron 
beam. The obtained yp center of mass energies are up to ten times the values reached by 
fixed target experiments. The integrated luminosity achieved in the first years of HERA 
running has produced a great yield of physics output. With increasing luminosity the 
improvement of some of these measurements relies on a quite precise knowledge of the 
responses of the various detectors. The HERA results have also triggered a huge amount 
of theoretical activity. 

The held covered by photoproduction at HERA has several connections with physics at 
e+e” and hadron colliders. The photon structure is probed in a complementary way with 
respect to the deep inelastic ey scattering at colliders. Moreover most of the photo¬ 
production processes are seen to resemble hadron-hadron processes and can be described 
in a similar framework. It is non-trivial to say that we are testing our understanding 
of the high energy scattering with different beam particles and the universality of the 
underlying picture. 

We review here the photoproduction results from HERA. In section ^ an introduction 
is given on the concept of the photon as a hadron, the measurements of its hadronic 
structure in the ey deep inelastic scattering, the definition of parton distributions for the 
photon and the most used sets of parametrizations. The basic characteristics of HERA, its 
kinematics and the experimental procedure to select photoproduction events are described 
in section ^ and a brief description of the ZEUS and HI detectors is given in section |^. 
The physics topics are then grouped in the following way: soft physics is dealt with in 
section |^, where the measurement of the total yp cross section and its decomposition 
in the different diffractive and non diffractive processes are described and in section |^, 
where the elastic vector meson production is described; hard photoproduction is then the 
subject of section ^ where the emphasis is on both tests of the QCD dynamics and on 
the determination of the photon and the proton structure from measurements of inclusive 
particle, jets, open charm and inelastic production; section || is devoted to hard 
diffraction, where both photon dissociation with the appearance of jets or heavy flavours 
and colour singlet exchange are described. Section concludes with some remarks about 
the future prospects. 


2 Resolved photon and photon structure functions 


In the Standard Model the photon is the gauge particle exchanged in electromagnetic 
interactions. It is elementary and its fundamental interactions occur with the charged 
leptons and quarks. On the other hand hadrons are composite objects, bound states 
of quarks and gluons, and have a size of the order of 1 fm cm). Despite this 

classification the photon can also behave like a hadron. The dual nature of the photon 
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was first observed in fixed target photoproduction experiments using real photon beams 
m]. These experiments have covered center of mass energies from 1 to about 20 GeV. The 
hadronic behaviour of the photon has been interpreted in the Vector Dominance Model 
(VDM) [^. It describes the photon-nucleon interaction as a two stage phenomenon: 
firstly the photon fluctuates into a vector meson state with the same quantum numbers 
= 1 0), then a normal hadronic interaction takes place between the nucleon 

and the vector meson. The total cross section is written as: 

_ Trrv 

a(^N^X)= Y, —criVN^X) ( 1 ) 

V=p,uj,(j) '^v 

where 'yv is the photon-meson coupling, which can be measured from electromagnetic 
decay rates like V e~^e~. Photoproduction cross sections are suppressed by the QED 
coupling a in comparison to purely hadronic ones, but the characteristics of the final state 
are found to be the same as e.g. vriV interactions. 

The bulk of hadronic interactions are soft processes involving low momentum transfers. 
However at high enough energies hard interactions can occur, their signature being the 
production of high px particles or jets. In Leading Order of QCD (LO) the strong coupling 
constant has the expression: 

, » ( 2 ) 

(33 - 2nf) In (f,) 

with rif the number of active flavours. Whenever the characteristic energy scale p of the 
interaction, like the px of the produced jet, is much higher than the QCD scale A (?« 0.2 
GeV), as is small enough to apply perturbation theory. 

In hard jN interactions the photon can behave as a pointlike particle in the so-called 
direct photon processes or it can act as a source of partons, which then scatter against 
partons in the nucleon, in the resolved photon processes. Resolved interactions stem from 
the photon fluctuation to a qq state (or a more complex partonic state). Depending on 
the relative transverse momentum between q and g, this may be a bound state like in the 
VDM {hadronic component) or a perturbative calculable state {anomalous component). 
Both are embedded in the definition of the photon structure functions. 

The apparent hadronic structure of the photon can be studied in ey deep inelastic 
scattering at e^e“ colliders. This proceeds through the exchange of two photons, one 
highly virtual (probe) and the other almost real (target), as in Fig. 2. Experimentally 
one electron is detected at a large angle with respect to its initial direction, and its energy 
Etag and polar angle Otag are measured; the other electron is required to scatter at small 
angle, generally by an antitagging condition, and it is thus lost into the beam pipe (single 
tag). Without tagging the electron scattered at low angle, the energy and small virtuality 
of the target photon are undetermined. It is then necessary to measure the hadronic 
invariant mass Whad- The following kinematic variables are dehned, with reference to 
Fig. a {k is the four-momentum of the incoming electron): 
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( 3 ) 


The reaction e’''e“ —> e^e~hadrons can be regarded as e 7 —>■ e hadrons factorizing 
the flux of quasi-real photons, i.e. the lower vertex in Fig. |I[ This is achieved in the 
Weiszacker-Williams approximation or better in the Equivalent Photon Approximation 
Naming z the energy fraction carried by the quasi-real photon, the differential cross 
section versus integrated over the photon virtuality between the limits 

p 2 is- 

^(ee ^ eeX) = f^/e{z) ^(ey ^ eX) (4) 

with the photon flux given by: 
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Here Pi,„ is the kinematic limit: 


P^ . 


mlz'^ 

1-z 


( 6 ) 


while P^ax corresponds to the maximum scattering angle dmax for fhe electron not to be 
detected (typically dmax ~ 30 mrad): 


p2 


(1 - ^)ELm^ 


2 

max 


(7) 


The level of accuracy of the Equivalent Photon Approximation goes from few per cent to 
about 10 %, depending on the P^ range of integration. 

The cross section for the deep inelastic ey scattering is written in complete analogy 
to the electron-nucleon DIS, introducing the two photon structure functions P^ and F'l'. 


da^ej —>■ eX) 
dEiag dcosOfag 


Ana^Etag 

Q^y 


{[l + {l-yf]E^{x,Q^) 


y^E2{x,Q^)} 


( 8 ) 


In the usual experimental conditions y‘^ is rather small and P^ can be neglected. The cross 
section measurements thus correspond to determinations of the E 2 structure function. 

The most striking feature of P^ is that in the region of asymptotically high it is 
completely calculable by perturbative QCD (at large x). This has to be contrasted with 
the nucleon structure functions, for which QCD predicts for example the evolution 
but not the absolute normalization at any value. 
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The lowest order contribution to is given by the Quark-Parton Model diagram 
shown in Fig. || (plus its crossed one), which is purely electromagnetic. It gives P: 
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where the sum runs over all quark and antiquark flavours q. The charge (in units of the 
proton charge) and mass are given by Cq and niq. This expression depends on the quark 
masses and roughly describes experimental data using constituent masses of few hundred 
MeV. However it shows some important qualitative features. Comparing it to the known 
behaviour of the nucleon structure function two major differences appear: 


• There is no scale invariance holding for even in the Parton Model. For nucleons 
the scaling violations result from gluon radiation and have opposite signs in different 
X ranges: with increasing F^ increases at low x while it decreases at high x. 
Instead F^ grows at all the x values with increasing Q^. 

• F 2 is large at high x, while counting rules predict for the nucleon a vanishing F^ 
for a; —>■ 1 . 


The different behaviour of photon and nucleon structure functions originates from the 
pointlike coupling of the photon to quark-antiquark pairs 7 —> qq: this coupling is what 
endows the photon with a hadronic structure. Perturbative QCD can be applied on top 
of the Parton Model and resummations have been performed for F 2 both in leading 
and next-to-leading (NLO) order []^ of QCD for large x and asymptotically high Q^. The 
result reads as: 


F^{x,Q‘^) = a 


a(x) 

as(Q^) 


+ b(x) 


( 10 ) 


with a and b calculable functions of x but diverging for x —> 0. The first term is the LO 
result, the second its NLO correction. The absolute normalization of F^ would then be a 
direct measurement of asiQ"^), or equivalently of the QCD scale A. This result accounts for 
only pointlike couplings and is valid at high where the incalculable hadronic component 
of the photon is not important. Unfortunately this regime is not experimentally reachable 
at present and foreseeable future accelerators. Non-perturbative contributions due to the 
formation of bound states after the 'y —>■ qq splitting cannot be neglected. 


The successful approach to F 2 , firstly suggested by Gliick and Reya |^, requires the 
introduction of boundary conditions, which have to take inputs from experimental data at 
some arbitrary scale Qq already in the domain of perturbative QCD. At the scale Qq both 
perturbative and non-perturbative contributions are buried together. From the boundary 
conditions the evolution with Q^ is then described by perturbative QCD, most practically 
with generalized DGLAP equations p. At the starting scale Ql parton densities are 
defined with a quark-antiquark component g'^(x, Qg) = (p{x, Ql) and a gluon component 
g'^{x,Ql). The evolution equations for the photon parton densities differ from those of 
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the nucleon due to the anomalous coupling 7 —> qq. Its splitting function is given by the 
coefficient of the leading logarithmic term in (|^): 


ki{x) 
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The anomalous coupling introduces an inhomogeneity in the evolution equations, which 
are written in the form: 
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Here t = In , q] and g'^ are the quark (antiquark) and gluon density in the photon and 
Pij {i,j = q,g) the ordinary splitting functions. Equations (JT^) hold in this form at both 
LO 1^ and NLO [0, |T2|, ^], provided the right one or two loop splitting functions 
and coupling constant ag are used. Recently a full NLO calculation of the heavy quark 
contribution to the photon structure functions has been carried out [jT4|. 


A peculiar characteristic of the photon structure is that there is no momentum sum 
rule constraining the photon parton distributions. In contrast in the nucleon case the 
early DIS data showed that gluons should carry about half of the nucleon momentum, 
thus fixing the size of g^{x,QQ). The reason is that in the quantum decomposition of 
the real photon the hadronic states and thus the parton densities are of order aem, while 
the bare photon state is of zeroth order. Thus any variation, even large, in the parton 
densities can be reabsorbed in a renormalization of the bare photon wavefunction ||I5| , p! 6 | . 


Quark and antiquark distributions are however constrained by the measurements of 


F 2 , since in LO: 


nf 


{x, QQ = ^ xq] (x, QQ + xq] (x, 


2 = 1 


(13) 


The first data on F 2 came from experiments at the e’''e“ colliders PETRA and PEP 
(for a review see |^). Most of the available photon parton distributions have been 
parametrized from those data. Later results have come from the TOPAZ |^] and AMY 
|T^ experiments at TRISTAN and from the OPAL [2^, ^ and DELPHI [^| experiments 
at LEP. They are shown in Fig. compared to LO parametrizations obtained from pre¬ 
vious data. Most of the data points suffer from quite large errors. One of the limitations 
are the low statistics in comparison with similar measurements of the nucleon structure 
functions. The reason is that the two-photon exchange implies a cross-section of order 
instead of as for lepton-nucleon DIS. Another important contribution to the ex¬ 

perimental errors is of a systematic nature. As we said, the measurement of x requires 
that of Whad- Due to the limited acceptance of the detectors, particularly around the 
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beam pipe, one is not sensitive to the whole hadronic hnal state. To correct for detector 
effects a model of the hadronic hnal state is necessary and this is usually done with the 
help of Monte Carlo generators. Unfortunately the unfolding of the true Whad from the 
measured one is strongly dependent on the different models, especially at low x. The 
expected accuracy of the measurements which can be done at LEP2 in an enlarged kine¬ 
matic domain is certainly very desirable . The knowledge of the photon structure is of 


utmost importance at LEP2 for the research of new physics, given that the 77 processes 
are an important source of background. 

The logarithmic increase of F 2 with predicted by (^, is clearly established from 
experimental data in Fig. where F 2 is averaged in the large-x region 0.3 < x < 0.8. 
This behaviour is a proof of the anomalous component of the photon. To account for the 
normalization a non-perturbative VDM component is also needed, like in the FKP model 


2 ^ shown by the curves in the figure. 


Most of the current sets of parton distributions of the photon follow the approach of 
evolution equations. A deeper description can be found in |^. The Drees and Grassie 
parametrization (DG) was obtained by LO evolution equations starting from input 
parton distributions ed, Ql = 1 GeV^, to fit only one data point at that time available 
at = 5.9 GeV^. Its input consists of three quark flavours with densities assumed 
proportional to the squared quark charges and flavour thresholds for charm and beauty 
crossed through the evolution. The input gluon is generated radiatively and does not 
enter in the £t. 


A similar strategy was used in the parametrizations by Levy, Abramowicz and Gharchula 
(LAG) [^, which were based on the world data set available in 1991. However they in¬ 
cluded in the £t the gluon distribution at the starting scale, demonstrating that data on 
F 2 do not constrain it signihcantly. This is connected with the lack of a momentum sum 
rule for the parton distributions in the photon. Among the LAG sets, LAGl and LAG 2 
have a very soft gluon distribution, rising steeply at low x, while LAGS has a very hard 
distribution, with a maximum at x ~ 0.9. The LAG parametrizations are LO and use 
four massless quark flavours. LAGS has already been excluded by data on jet production 
in 7 p and 77 scattering. 

The WHIT parametrizations (after Watanabe, Hagiwara, Izubuchi and Tanaka) 
follow a similar approach as LAG, but treating correctly the charm threshold. This is 
more important than in the nucleon case, because the photon can easily develop a charm 
content through the anomalous splitting 7 —>■ cc. Six sets have been provided, each one 
with a different assumption on the shape of the gluon distribution. The total momentum 
carried by gluons is however fixed at about half the momentum carried by quarks (at 
Qq = 4 GeV^) in WHIT 1,2,3, while gluons and quarks carry about the same momentum 
fraction in WHIT 4,5,6. 

The Gordon-Storrow parametrizations (GS) exist at both LO and NLO. Their hrst 
version starts from a scale Qq = 5.3 GeV^ with a VDM ansatz modihed by the addition 
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of a pointlike term: 


Aircy 

{q\9^){Ql) = k^{q^ 

J O 




(14) 


Here the structure of the vector meson is assumed equal to the pion structure. The free 
parameters are the momentum fractions carried by gluons and sea quarks in the pion, 
the constant k and the light quark masses from the QPM term. In the GS2 distribution 
comes entirely from the first term, while in GSl the second term also contributes 
through radiation from the pointlike quark component. Recently these distributions have 
been updated by lowering their starting scale to 3 GeV^, including all available data 
on F 2 and constraining the gluon from jet production data at TRISTAN. 


The GRV parametrizations (after Gliick, Reya and Vogt) exist both in LO and 
NLO and are built at a very low starting scale: Qg =0.25 GeV^ (LO), 0.3 GeV^ (NLO). 
Such a low value is considered the minimum for the applicability of perturbative QGD. 
The authors apply the idea of dynamically generated parton distributions, which they 
had previously used to describe the nucleon and pion structure [^. The ansatz for 


the nucleon was originally to have only valence quarks at the initial scale. It was 
then modihed introducing valence-Iike distributions of gluons and sea quarks. For the 
pion no sea quarks are needed at the starting scale Ql but a valence-like gluon density, 
proportional to the valence quark density. The GRV approach has attained considerable 
success with HERA data on the proton F 2 structure function [^, although it is not 
obvious that perturbative QGD implemented in DGLAP equations works already at < 
1 GeV^. The GRV parametrizations for the photon start from a VDM input: 


{q\9^){Ql) 


Atiq. 

J O 




(15) 


where only a normalization factor k {1 < k < 2) remains free, accounting for p/uj/cf) 
mixing. In this way the gluon distribution at the scale Qg is constrained to be equal to 
that of the pion. 


The AFG parametrization (after Aurenche, Fontannaz and Guillet) is built in 
NLO and, like GRV, starts from a low scale Qg = 0.25 GeV^ assuming a VDM ansatz. 
However a coherent sum of p, a;, 0 is supposed. Pion distributions are taken from . 


The authors choose the MS regularization scheme pointing out that the DIS^ scheme 
adopted by GRV is more stable but includes process dependent terms in the parton 
distributions. 


Schuler-Sjostrand parametrizations are quite recent and follow a similar approach 
as GRV and AFG. They exist only in LO. Different sets are provided, including parton 
distributions for the virtual photon. 


Summarizing, the VDM ansatz allows one to constrain the gluon distribution in the 
photon, otherwise free due to the lack of a momentum sum rule. Unfortunately with 
typical starting scales Qg >1 GeV^ a pure VDM input is not able to reproduce data at 
higher Q^. The various parametrizations solve the problem by adopting one of these two 
approaches: 
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• starting from a very low scale Qq with a pure VDM input, in a region where per¬ 
turbative QCD may be doubted; 

• keeping the scale at Qq >1 GeV^ and using a more complex input with the addition 
of a pointlike component given by the Quark-Parton Model or simply guessed. 


A comparison between the quark distributions from different sets is shown in Fig. ^ There 
is reasonable agreement between them in the region 0.05 < x < 0.8, where data exist. 
Instead large differences appear at high x, for example comparing GS and GRV. The 
striking difference between the quark distribution in LO and NLO is characteristic of the 
MS scheme, and this has been discussed in [^. It does not affect the hnal predictions on 
observable quantities like F^. It can however be avoided using the DIS^ scheme which is 
adopted by the GRV parametrizations . 


Large differences between the parametrizations also occur at low x. This is clearly 
visible in Fig. p, where different LO gluon distributions are compared. The LAG sets 
are much steeper at low x than the others, demonstrating that the ey DIS data do not 
constrain the gluonic component of the photon. At NLO the similarity of the curves is 
due to the similarity of the models, it is not borne out by data. 


The hrst experimental constraints on the gluon density in the photon have come from 
experiments at TRISTAN (a/s = 58 GeV) from the study of 77 interactions, with 
both photons almost real. Here jet production occurs in LO through the Quark-Parton 
Model process as in Fig. [^a, or through single- or double-resolved processes as in Fig. |^/b- 
c. The three kinds of processes are all of the same order in the coupling constants since 
the photon parton distributions are 0(aem/«s)- The first experimental evidence for a non 
zero gluon content came from AMY [^]. Single jet and dijet inclusive cross sections have 
then been measured by both TOPAZ and AMY |^. The results of TOPAZ compared 
to theoretical predictions are shown in Fig. H. The need for the hadronic component of 
the photon is evident when comparing data to the expected cross section for direct photon 
only. Moreover it is clear that there is a need for a gluon component in the photon. These 
data also rule out the LAGS parametrization which would predict a far too high cross 
section. 


3 HERA as a collider 


The HERA storage ring at DESY (Hamburg, Germany) is the hrst electron-proton collider 
ever built. Two different magnetic systems, respectively superconducting for protons and 
conventional for electrons, drive the beams in separate rings along a circumference of 6.3 
Km. The design energies are = 30 GeV, Ep = 820 GeV, giving a center of mass energy 
a/s = 314 GeV (currently Eg = 27.5 GeV, so that ^/s = 300 GeV). This corresponds to 
an electron beam energy of 52 TeV against a hxed target. The beams cross each other 
at zero angle in two interaction regions, occupied by the experiments HI and ZEUS. The 







bunch crossing occurs every 96 ns. The designed number of bunches in each beam is 210. 
To check beam related backgrounds about 10 bunches are left unpaired in each beam. 

The HERA kinematics is represented in Fig. p. The invariant kinematic variables of 
the semi-inclusive reaction ep —*• eX are: 


Q" = 

—q^ = —(k — k'Y 

(16) 

X = 

2P-q 

(17) 

y = 

P ■ q 

P-k 

(18) 

= 

{P + qf 

(19) 


Two of them can be chosen as independent variables. is the virtuality of the exchanged 
photon. The photoproduction regime is dehned by low Q^, such that the photon is almost 
real. In this case we may regard HERA as a photon-proton collider, where the photons 
are emitted almost collinearly with the electron beam. The photon spectrum can be again 
evaluated with the Equivalent Photon Approximation, factorizing the vertex e —> ey in 
the flux factor of formula with the kinematic limits given by (^) and (|^), in which the 
variable is substituted by y. Due to the asymmetric energies of the HERA beams the 
7 P center of mass is boosted typically by two units of pseudorapidity y, dehned as: 


j P + Pz 
In - 
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—In 



( 20 ) 


This means that a particle going at 90° in the yp center of mass appears at about 15° in 
the HERA laboratory frame. In the following sections we will refer to forward or backward 
(rear) directions intending respectively the proton and the photon directions. 

The photoproduction regime is experimentally dehned in two ways: 


• tagged photoproduction, when the scattered electron is detected at a very small 
angle. This can be done in both ZEUS and HI by dedicated calorimeters placed 
near the beam line at several meters from the interaction point. The acceptance 
extends in this case from 10“® GeV^ to 0.02 GeV^; 

• untagged photoproduction, dehned by the antitagging condition that the scattered 
electron is not detected in the main calorimeter. This dehnition has been used 
by both the experiments. It corresponds to < 4 GeV^, with median about 
10-^ GeVT 


In contrast the region of > 4 GeV^ is the traditional Deep Inelastic Scattering 
(DIS) region, where the process in Fig. ^is understood as a pointlike interaction between 
the electron and a quark constituent of the proton. In this regime the structure of the 
proton can be studied by exploiting the almost pointlike nature of the exchanged virtual 
photon. This region is experimentally dehned by tagging the scattered electron in the 
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main calorimeter. The x variable represents in this context the momentum fraction of the 
struck quark relative to its parent proton. 

The meaning of y is that of scaled photon energy: in the proton rest frame 0 < 
y = E^/Ee < 1 and W is the photon-proton center of mass energy or equivalently the 
invariant mass of the hadronic final state. In the photoproduction regime these are related 
by IT ~ a/^- III comparison with previous photoproduction experiments, carried out 
with fixed targets, HERA enlarges by one order of magnitude the range of center of mass 
energies. Depending on the specihc process the HERA experiments can span the whole 
range of W between 30 GeV and 300 GeV running at the nominal beam energies. Lower 
values of W could be accessed by reducing the energy of the proton beam. 

In the tagged mode y (and hence W) is easily determined by measuring the energy of 
the scattered electron E' in the low-angle electron tagger by: 

2/ = 1 - f (21) 


In the untagged mode y can be determined from the measurement of the hadronic hnal 
state by the Jacquet-Blondel variable pi|, dehned as: 


_ E.^^(l-cosg,) 
{E - p^)tot 2Ee 


( 22 ) 


where the sum runs over all the calorimeter cells, Ei being the energy contained in the i-th 
cell at polar angle 9i with respect to the proton direction, yjs gives a systematically lower 
value than the true ?/, mostly due to particles escaping detection inside the beam pipe in 
the photon direction or inactive material in front of the calorimeter, and is corrected by 
Monte Garlo detector simulation. 


The hnal states produced in photoproduction events can be very diverse, from hadron- 
like minimum bias yp interactions to (quasi)elastic vector meson production, to the pro¬ 
duction of large pr jets and heavy havours. Detection of all those processes demands 
multipurpose detectors such as ZEUS and HI with many capabilities and good overall 
performances. 


4 The experiments ZEUS and HI 

Both ZEUS and HI have multipurpose detectors. Here only a brief description is given, 
more details can be found in |43[| . 

A longitudinal section of the central part of the ZEUS detector is shown in Fig. 0. 
Gharged particles are measured by the inner tracking detectors: the vertex detector 
(VXD), the central tracking detector (GTD) and forward and rear tracking detectors 
(FDET, RTD). They operate in a magnetic held of 1.43 T produced by a thin super¬ 
conducting solenoid. The VXD and the GTD are concentric cylindrical drift chambers 
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covering the angular region 15° < 6 < 164° (where 6 is the polar angle with respect to the 
proton direction). The measured resolution for tracks going through all the CTD layers 
is: <Jpj,/pT = 0.005 p'r(GeV)0 0.016. 

Surrounding the coil is the high resolution uranium-scintillator calorimeter (CAL), 
mechanically divided into three parts, the forward (FCAL), barrel (BCAL) and rear 
(RCAL) calorimeter. The calorimeter is compensating, the thickness of the uranium and 
scintillator tiles were optimized to achieve the best possible energy resolution for hadrons. 
It covers the angular region: 2.6° < 6 < 176.1°, corresponding to 99.7% of the solid angle, 
with holes of 20 x 20 cm^ in the center of the forward and rear calorimeters to accomodate 
the HERA beam pipe. Each part of the calorimeter is subdivided into towers of transverse 
size 20 X 20 cm^ which are segmented longitudinally into one electromagnetic (EMC) and 
one or two hadronic (HAC) sections. In the electromagnetic section the readout cells have 
finer granularity of 5 x 20 cm^ (in FCAL and BCAL) or 10 x 20 cm^ (in RCAL). Each cell 
is read out by two photomultipliers. Under test beam conditions the energy resolution is 
ge/E = 0.18/v^ for electrons and ge/E = Q.^b/\fE for hadrons, with E in GeV. The 
time resolution is better than 1 ns for energy deposits greater than 4.5 GeV. 

Outside the uranium calorimeter a moderate resolution calorimeter (BAG), made of 
iron layers interleaved with proportional tubes, measures the tails of high energy jets 
and acts as a hlter for muons. The muon detectors are also divided into three sections, 
covering the forward, the barrel and the rear regions. In the barrel and rear regions large 
chambers made of limited streamer tubes measure the position and the direction of muons 
in front and behind the iron yoke (BMUON, RMUON). The yoke is magnetized to 1.6 T 
to allow the momentum measurement for penetrating muons. In the forward direction a 
spectrometer (FMUON), made of two iron toroids, drift chambers and planes of limited 
streamer tubes, identihes muons and measures their momenta. 

Charged particles scattered at very small forward angles and carrying a substantial 
fraction of the proton beam momentum are detected in the leading proton spectrometer 
(LPS). This is a set of six stations placed along the beam line at distances from 20 to 90 
meters from the interaction point, using planes of silicon microstrip detectors very close 
to the proton beam. The track deflection induced by the magnets in the proton beam 
line is used to measure the momentum. The LPS measures a leading proton for fractional 
momentum xl > 0.4 and pr < 1 GeV. 

The luminosity is measured from the rate of the Bethe-Heitler process ep eyp, whose 
cross section is large and calculable to high precision. Two lead-scintillator calorimeters 
placed downstream the electron beam pipe detect respectively the outgoing electron and 
the photon. The electron calorimeter lies at 35 m from the interaction point and accepts 
electrons with energy between 7 and 20 GeV scattered at angles up to about 5 mrad. The 
photon calorimeter is placed at 107 m from the interaction point and accepts photons 
produced at angles below 0.5 mrad. Bremsstrahlung events are tagged by the coincidence 
of the two calorimeters where the energies add up to the energy of the electron beam. 
The luminosity is calculated from the photon tagger alone. 
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Figure p!T| shows a schematic view of the HI detector. Charged particles are tracked by 
the central and forward tracking systems, both consisting of drift and multiwire propor¬ 
tional chambers. The momentum resolution achieved is <Jpj,/pT ~ 0.009■p'r(GeV) 00.015 
and the angular coverage is the same as for the ZEUS detector. 

The tracking region is surrounded by a hne grained liquid argon calorimeter (LAr) with 
an electromagnetic section made of lead absorber plates and an outer hadronic section 
with steel absorber plates. Its angular coverage is 4° < 6^ < 153°. The energy resolution 
under test beam conditions is (Je/E = 0.12/a/E for electrons and aE/E = O.hO/A/E for 
hadrons, with E in GeV. The calorimeter is not compensating and offline corrections are 
applied according to the shower shape and energy. Its segmentation varies as a function 
of the polar angle and is hner in the proton direction where, due to the asymmetry of the 
HERA beam energies, the particle density is higher. In this region the transverse cell size 
is as small as 3 x 3 cm^. The backward region 155° < 6^ < 176° is covered by another 
calorimeter, which was primarily intended for the detection of the scattered electron in 
deep inelastic scattering. This consisted of a lead-scintillator sandwich with depth of 
22.5Xo and energy resolution (JeIE = O.IO/a/E. From the 1995 data taking it has been 
replaced by a scintillating hber calorimeter with photomultiplier readout. 

The calorimeters are contained inside a superconducting solenoid providing a uniform 
magnetic held of 1.15 T for the inner tracking detectors. Outside the coil there is the 
iron yoke, which is instrumented with LST planes for the muon identihcation. Further 
muon chambers in front and behind the iron complete the muon detectors in the central 
region, while in the forward region muons are measured in a toroidal spectrometer similar 
to the one in ZEUS. The luminosity is measured in a similar way as in ZEUS, from the 
coincidence rate of two TlCl/TlBr crystal calorimeters. 

HERA operation began in summer 1992 with electron beam energy of 26.7 GeV. Since 
1994 positrons are run at energy of 27.5 GeV. The luminosity has continuously increased 
from year to year. In Fig. the integrated luminosity delivered by HERA is shown as 
a function of time for the period 1992-1995. We review here the hnal results available at 
present, which have been obtained analyzing data taken up to the 1994 run. 


5 Total 7 p cross section and its decomposition. 


The high energy behaviour of total hadronic cross sections is related to fundamental 
properties of particle interactions, pp colliders have provided data up to a center of mass 
energy ^/s = 1.8 TeV, demonstrating the logarithmic rise of the cross section ^ ^ . 
Asymptotically this rise cannot violate the Froissart-Martin bound [^, Q, which states 
that total cross sections cannot grow faster than ln‘^{s), as a consequence of unitarity. At 
HERA the yp cross section has been measured up to a center of mass energy of 200 GeV, 
an order of magnitude higher than previous photoproduction measurements from hxed 
target experiments . 
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The most popular approach to hadronic cross sections employs Regge theory ideas 


They have been able to £t all the available cross section data (for ^/s > 5 GeV) using a 
parametrization inspired by Regge theory, of the form: 

o'tot = + Rs ^ (23) 


M. In this framework the simplest description is due to Donnachie and Landshoff K 


The hrst term arises from pomeron exchange and the second from reggeon exchange (p, 
uj, f, a). The exponents e and 6 are effective powers, summing effects of multiple pomeron 
or reggeon exchanges. Such effects are supposed to be small, therefore e and 6 are related 
to the intercepts of the pomeron and reggeon trajectories: 


— Q;jp(0) -|- cxjp ■ t 

(24) 

a-Tzit) = a-jiiO) + 

(25) 

e = ajp(O) - 1 = 0.0808 

(26) 

5 = 1- an{^) = 0.4525 

(27) 


that is: 


The pomeron has the quantum numbers of the vacuum, so that its coupling to a particle 
a and its antiparticle a are equal. Thus the coefficient X is set equal for a{ab) and a{db). 
This is a way of enforcing the Pomeranchuk theorem in the fit. The universal rise 


of the cross sections with energy is interpreted by Donnachie and Landshoff as sign of 
the universality of something that is exchanged, the pomeron, which is sometimes also 
referred to as the “soft” pomeron. This simple parametrization becomes inconsistent at 
very high energies, violating the Froissart-Martin bound. However this takes place only 
at ^/s > 10^^ GeV, very far from the presently conceivable energies. 


On the other hand several models have been based on perturbative QGD, attributing 
the rise of the cross section to the onset of “minijets” |^. In this context the total 
inelastic yp cross section is assumed to have a energy independent part, coming from soft 
non-perturbative interactions, and a QGD part rising with energy, due to the increase of 
gluons in the proton and in the photon at low partonic energy fractions. Predictions based 
on this approach depend on the parton distributions of the proton and the photon, the 
minimum transverse momentum down to which the perturbative calculation is assumed 
to be valid (p™""), and the probability for a photon to go into a hadronic state. These 
models use the eikonal formalism to ensure unitarity and in general have a large 
number of unknown parameters. 


The total yp cross section is the sum of quite different processes, represented in Fig. |T^. 
Let us first divide diffractive (|T^a-d) from non-diffractive processes (P!3|/e). In diffractive 
processes there is no colour exchange between the proton and the hadronic state of the 
photon. In the figure they are represented to proceed through pomeron exchange. As¬ 
suming the validity of the VDM for the transition j —>■ V {V = p, u, (p) the following 
four classes are identified: 
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• elastic vector meson production ■yp ^ Vp (Fig. p!3|/a). The truly elastic Compton 
scattering yp —>• jp has a cross section about two orders of magnitude smaller and 
can be neglected. 

• photon diffractive dissociation jp —>■ Xp (Fig. ^/b). In this case the photon gives 
rise to a low mass hadronic system X with the same quantum numbers, while the 
proton stays almost untouched and is the leading particle in the final state. 

• proton diffractive dissociation yp —>■ VY (Fig. m/^)- The photon converts to a 
vector meson and the proton breaks up into a low mass hadronic system Y. 

• double diffractive dissociation yp ^ XY (Fig. [T^/d). Here both the photon and 
the proton dissociate to low mass hadronic systems. 


In all cases the transverse momentum exchanged is generally low and the topology of the 
final state is characterized by large rapidity intervals (or gaps) without particles. This is 
expected if there is no colour flow between the interacting particles and is exploited as an 
experimental signature. 

Non-diffractive processes (Fig. are mostly soft ones, similarly to minimum bias 

events at hadron colliders. In this case colour exchange occurs and thus final state particles 
are distributed throughout the available phase space without signihcant gaps. 

It should be clear that the different processes listed above are qualitatively quite dif¬ 
ferent from one another. Experimentally, to trigger and select inclusive photoproduction 
events at HERA one needs a coincidence between tagging of the scattered electron in the 
low-angle electron taggers and activity in the main apparatus, at least a charged track or 
a minimum energy deposition. The detection efficiency for each of the photoproduction 
processes is quite different, thus complicating the measurement of the total cross section. 
Introducing the partial cross sections cxj (i=l,5) for the five considered classes of reactions, 
with cT(yp —Y) = cTj, the observed number of events Nev is given by: 

5 

Nev = ^ *f>7/e Ae 6^(7* (28) 

i=l 


with C the integrated luminosity, ^-^/e fhe photon flux, Ae the efficiency to tag the low 
angle electron and e* the efficiency of the main detector for the Fth process. The photon 
flux is obtained by integrating expression (|^) in the kinematic region corresponding to 
the low-angle electron tagger. The e* are calculated by Monte Carlo simulations of each 
different process. To reduce the systematic error it is important to determine separately 
each of the a,, otherwise one would have to rely on model dependent assumptions in the 
Monte Carlo, such as the relative magnitude of each contribution and its characteristics. 


The HI experiment has measured the total yp cross section by short dedicated 
runs with special triggers in 1994. They took 24 nb~^ data with nominal vertex position 
and the same amount of data with vertex displaced by about 70 cm towards the forward 
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part of the detector. In fact the acceptance for all the diffractive processes increases if the 
position of the interaction vertex is shifted in this way. The modest integrated luminosity 
is already enough due to the large photoproduction cross sections. The partial cross 
sections for the diffractive reactions are determined from hts to selected distributions, 
related to the rapidity gaps in the pattern of the hnal state particles. This is done on 
the shifted vertex data. The double dissociation cross section {(Tdd) is assumed to be 
in the reasonable range 0 — 40 fih. The non-diffractive cross section is then determined 
from nominal vertex data, giving as input in ( p8[) the diffractive partial cross sections 
previously determined. The measured cross sections are summarized in table 1. The 
elastic cross section aEL is almost independent of the assumed value of (Tod- Instead the 
single photon dissociation {(Tgd) and the single proton dissociation (apo) cross sections 
depend markedly on a^D, hut for any chosen value of apD the photon dissociation cross 
section is found to be larger than the proton dissociation one. The dominant systematic 
error (~ 5 %) in the result is on the acceptance of the electron tagger. This depends 
strongly on the electron beam optics, particularly on the offset and the tilt of the beam 
in the horizontal plane. 

The total cross section as measured by HI at an average c.m.s. energy W ~ 200 GeV 
is shown in Fig. together with previous data at lower energies and a measure- 

0 . 


ment from ZEUS 


The curves are predictions of Donnachie-Landshoff (DL) 


and Abramowicz-Levin-Levy-Maor (ALLM) 0 parametrizations, both of which did not 
include HERA data. The ALLM is based on both yp and 7 *p data at fixed target exper¬ 
iments. It allows a unified representation of photon-proton interactions with from 0 
up to the DIS regime. A DL fit taking into account the CDF measurement gives the 
exponent e = 0.11 in The HERA data are well accomodated in the soft pomeron 

phenomenology expressed by these parametrizations, thus supporting the VDM and the 
universality of the rise in the cross sections at high energy. 

Minijet models give predictions spread over a larger range due to uncertainties in 
the different parameters involved [Q. The present HERA results can be brought into 
accordance with calculations having a not too small p™"' (> 2 GeV) and a moderately 
rising gluon density in the photon, but strong conclusions cannot be drawn. 

In table 2 the diffractive partial cross sections measured by HI are compared to pre¬ 
dictions of several theoretical models, based on different assumptions about the structure 
and dissociation of the photon and the proton. CKMT |5^ is a calculation using Regge 
theory plus absorptive corrections for all the different processes. The measured cross sec¬ 
tions are derived fixing apo = 15 /rb, in agreement with the models. The elastic cross 
section is correctly predicted by all the models. The photon dissociation cross section is 
measured to be three times larger than the proton dissociation one, in agreement with 
predictions of GKMT. Instead SaS and GLM predict app > aoD- 
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6 Elastic vector meson production 


The elastic reaction 7 p —Vp (V = p, a;, 0, J/^, ...) represented in Fig. |T^/a is 
equivalent to a hadronic elastic interaction Vp — Vp, assuming the validity of the VDM. 
The optical theorem relates the total Vp cross section to the imaginary part of the forward 
elastic scattering amplitude: 




Im T{s, t = 0) 


(29) 


where t is the square of the momentum transfer. Soft pomeron phenomenology thus gives 
straightforward predictions on the elastic cross sections The elastic cross section is 


predicted to rise slowly with the c.m.s. energy and to be dominated by forward scattering, 
with a sharp peak at t = 0 of roughly exponential behaviour: 


da 

d\t\ 




(30) 


The 6-slope is related to the radius of the hadronic interaction and it can be understood 
in the context of a simple geometrical model. In a hxed target conhguration, for small 
angle of scattering t ~ —{p0)‘^, where p is the incident beam momentum and: 


— = ~ 1 - h{pef 

=0 

In the classical optical diffraction of light from a circular aperture of radius R the intensity 
of light is given by: 

L^i- ^(kef ( 32 ) 

where k is the wave number of the light. Comparing these two expressions one has a 
relation between the slope parameter and the radius of the interaction: 




b = 


T 


(33) 


For a typical radius of the strong interaction R = \/m.„ (with the pion mass) the 
last expression gives b ~ 12.5 GeV“^. This is in the measured range of values for elastic 
scattering of hadrons at high energy. The t spectrum is predicted to exhibit a dip followed 
by a secondary diffraction maximum, as for light diffraction. This is in fact observed in 
high energy pp and pp elastic scattering |^. Regge theory predicts a weak dependence 


of the 6-slope on the center of mass energy, referred to as shrinkage, of the kind: 


6 = 6o -|- 2a']p log — 
■So 


(34) 


This effect actually allows the determination of the slope of the pomeron trajectory in 
(H): a'jp ~ 0.25 GeV^. 
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Fixed target experiments have studied the elastic vector meson photoproduction up to 
W ^ 20 GeV. Light (p, to, 0) vector meson data EM have been found consistent with 
soft pomeron phenomenology and the VDM. Elastic J/'ip production cross sections 
were overestimated by VDM by more than an order of magnitude, but other predictions 
like the energy dependence agreed with data within the errors. 


UJ 


HERA experiments have measured the elastic reaction 7 p 




73, 74 


Vp for V = p 1^, 1^ , 

These measurements are complemented by similar ones 


carried out in the DIS regime ||75|, [76| . The reconstruction of the vector meson decays: 



—7r’'"7r 

UJ — 

•4^ 

+ 

1 

o 

(f) - 


J/iJ^ 

e~^e~, p~^p 


has allowed the measurements at W from 40 to 100 GeV for light mesons, from 30 to 150 
GeV for J/'0- One of the hardest experimental tasks is to trigger on such elastic events. 
In fact both the outgoing electron and the proton are generally not tagged. Then only the 
decay products of the meson decay can be observed in the detector. The typical topology 
of these events is of two charged tracks in the photon hemisphere (plus an electromagnetic 
calorimeter cluster due to the 7r° in uj decays). The acceptance in W is determined by the 
detector geometry and trigger. 


In tagged photoproduction events the p production has also been measured without 
reconstruction of the decay at (W) ~ 200 GeV, from fits to distributions of the hadronic 
hnal state by comparison to Monte Garlo distributions for p production |5^, ^ . At these 
high W values the decay pions are boosted in the photon direction outside the acceptance 
of the tracking detectors and their reconstruction is carried out by the calorimeter. In this 
case the worse spatial and energy resolutions have prevented measurements of differential 
distributions. 


The kinematic variables W and t can be easily determined from the vector meson, 
reconstructed from its decay products, in the limit of ~ 0. In fact in this case 
{pz)'y — —E^ and: 

~ AE^Ep ~ 2(£' — pz)vEp (35) 

t - -{Pt)v (36) 

where the subscript V indicates the vector meson. The main background to the elastic 
reactions comes from events where the proton dissociates (fig. p!3|/c). This can be reduced 
by requiring a very low energy deposit in the forward part of the calorimeter {E < 
1 GeV). The residual background is estimated studying the effect of changing this cut in 
comparison to Monte Garlo simulations of the relevant diffractive process. In the hnal 
samples it is generally 10 — 30 % and is statistically subtracted. 

Proton dissociation can be completely avoided by detecting the leading proton scat¬ 
tered at a very small angle. This has been done by ZEUS using the leading proton 


17 




































spectrometer (LPS), which directly measures the momentum of the leading proton with 
high accuracy (with the conhguration the detector had in 1994 run: 0.4% on the longitu¬ 
dinal component, ~ 5 MeV on p^)- For truly elastic events the momentum of the proton 
is reduced by less than 0.2% and thus events are selected with 0.98 < xl < 1.02, where 
xl = v'IV is the ratio between outgoing and incoming proton momentum {xl can be 
measured greater than unity due to the experimental resolution). One is thus guaranteed 
that the proton scatters elastically and t is directly measured. Owing to the energy- 
momentum balance the can also be determined in this case. In fact in the elastic 
reaction the knowledge of the pt of the vector meson and of the scattered proton allows 
a determination of the pt of the untagged scattered electron, and from this the Q^\ 


Q^ = 


Pt 


y 


(37) 


where y ~ W'^/s. The measured da/dt for elastic production is shown in Fig. A fit 
of the form (pOD gives b = 9.8 ±0.8 (stat.) ±1.1 (syst.) GeV“^, with dominant systematic 


uncertainties on the LPS acceptance and the beam spread in the transverse momentum. 

HERA measurements of the slope for p elastic photoproduction are shown in Fig. p±|in 
comparison with hxed target results. They agree with a slow shrinkage of the diffractive 
peak as predicted by the soft pomeron formula (|3^). This is also true for the measurements 
of the other light vector mesons, summarized in table Recalling that b is related to the 
radius of the interaction, the values in table indicate a smaller radius for (pp as compared 
to pp interaction. 

The integrated elastic cross sections for p, cn, 0, J/V’ photoproduction measured at 


HERA as functions of W are shown in Fig. together with data from hxed target 
experiments. In the soft pomeron model [Q one expects: 


a. 


ip- 


.Vp(W) 




4e 


OC 


b(W) 


~ W' 


0.22 


(38) 


where e is given by (PB]) and b(W) is given by (|5^ ) with s = W^. This is in good agreement 
with data on light vector mesons. The same model predicts for the total photoproduction 
cross section a dependence ~ which also agrees with data, as shown. 


Another characteristic of the elastic vector meson production at hxed target energies 
is the so-called s-channel helicity conservation (SCHC). For real or quasi-real photons it 
means that the transverse polarization of the photon is transferred to the vector meson. 
The polarization of the vector meson can be experimentally determined from the angular 
distributions of the decay products, which depend on the spin-density matrix elements 
771. The analyses reveal in all cases that SCHC continues to hold at the HERA energies. 


This homogeneous situation changes dramatically when one turns to the elastic pro¬ 
duction of J/4’- Coming back to Fig. 0 the J/'ip elastic cross section is steeply rising 
with W and the Donnachie-Landshoff model fails, predicting a slow rise as for the light 


vector mesons. In Fig. [Tq/ a only the ZEUS data are shown. The ht with a function W° 
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gives 6 = 0.92 ± 0.14 ± 0.10. Thus even in this restricted range the soft pomeron model 
is excluded by data. Instead other two pomeron models are found to give a good descrip¬ 
tion of the total cross section [^, ^ (Fig. |T^b). They involve either a scale dependent 
pomeron intercept or scale dependent pomeron couplings and mass threshold effects. The 
data can also be reproduced by models based on perturbative QCD 1^, where the 


pomeron is considered as a system of two interacting gluons, or gluon ladder, as in Fig. |T^. 
In these models the cross section is proportional to the square of the gluon density in the 
proton. For a singular gluon density at small x, rising like: 


X 


-A 


as X 


0 


(39) 


as implied by HERA F 2 measurements [^, a steep rise is predicted with increasing 
W. HERA F 2 data are well described by the MRS(A’) parametrization [Q, which cor¬ 
responds to a power A ~ 0.2 in (P^D, and by the GRV parametrization |Q corresponding 
to A ~ 0.3 — 0.4. In the Ryskin model fFe effective energy scale q and momentum 
fraction x at which the gluon distribution is tested are given by: 


f = 


-1 


X = 


i 


PP2 


(40) 

(41) 


For 0, f ~ 0 the only hard scale is represented by the J/ijj mass and the appli¬ 
cability of perturbative QCD relies on it. From 00 one has: ~ 2.5 GeV^ and 

5 ■ 10“^ < X < 5 ■ 10“^. Due to the quadratic dependence on the gluon density the shape 
of the cross section as a function of W is very sensitive to different proton parton density 
parametrizations. However the normalization of the theoretical result suffers presently 
from quite large uncertainties. In Fig. ^/b the elastic J/'ip data are compared to the 
prediction of the Ryskin model with input gluon density given by MRS (A’). A good 
agreement is found, while the steeper behaviour of GRV has been reported to overshoot 
HERA data []75| . 


HI has determined the slope of the J/il) distribution by htting data with < 1 
GeV^ to an exponential in the range 30 < IE < 150 GeV, achieving h' = 4.0 ± 0.2 ± 
0.2 GeV“^. The error in the slope due to using p^ instead of t is evaluated from Monte 
Carlo to be —10%. No clear evidence for shrinkage is obtained dividing the sample in W 
bins, mainly due to the still limited statistics. A consistent result comes from the ZEUS 
analysis, which also corrects explicitly the p\ differential cross section to obtain da/d\t\, 
shown in Fig. Fitting over the range \t\ < 1 GeV^ with an exponential shape 

lo;g GeV“^, which is consistent with the slope values from fixed target 

ill- 


gives 6 = 4.6 ± 0.41q g 
experiments 


HI also measured the cross section for production with proton dissociation, ex¬ 
ploiting their forward detectors. This turns out to be of the same order of the elas¬ 
tic cross section, in a similar range of W. The p^ distribution gives a slope: b' = 
1.6 ± 0.3 ± 0.1 GeV“^, a factor of 2.5 smaller than for pure elastic scattering. 
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The study of the decay angular distributions of the J/ijj (both elastic and with proton 
dissociation) shows consistency with s-channel helicity conservation, i.e. the J/'0, like 
lighter mesons, is mainly transversely polarized. 


Diffractive production of 'ifj' has been observed by HI in the decay chain ip' —> 
JI'll) 7r+7r“, J/pj ^ The rate of its production compared to J/pj is found to be 

0.20 ± 0.09, in agreement with previous measurements at hxed target experiments. 


Elastic vector meson production has also been measured at high for 0 and 
J/p) 1^, Preliminary results exist for low and high-t photoproduction with proton 
dissociation ||^ . These data, compared to fixed target data, suggested a hard production 


mechanism when a high energy scale is present in the interaction, be it a high mass (as 
the J/p) mass) or high or t. However recent results from the E665 Collaboration 


appear at variance with the previous hxed target data. The experimental situation 
for light vector mesons is still to be clarihed. Thus vector meson production at HERA 
allows to study the transition from a non-perturbative description of the pomeron (the 
soft pomeron), to the hard perturbative pomeron. 


7 Hard photoproduction and QCD 


Hard photoproduction processes at HERA are of particular interest. The large cross 
sections and the range of high center of mass energies allow precise studies of the photon 
structure in a way complementary to the ey deep inelastic scattering. Moreover the QCD 
dynamics and the peculiar features of photon induced reactions can be tested. The hard 
scale is provided by the highest transverse momentum of the produced particles or jets in 
the final state. Application of perturbative QCD requires that this scale is at minimum a 
few GeV. As mentioned earlier, in leading order of QCD two classes of processes contribute 
to the production of high pj- particles or jets. The photon may interact directly with a 
parton in the proton, as in Fig. ^/a, or it may hrst fluctuate into a hadronic state as in 
Fig. |2T|/b and then act as source of partons, which scatter off partons in the proton. In 
the hrst case, known as direct photon, the whole photon momentum enters in the hard 
subprocess and the hnal state consists of two jets balancing each other in pt plus the 
proton remnant and the scattered electron, emerging with small pt- In the second case, 
known as resolved photon, only a fraction of the photon momentum participates in the 
hard process and the hnal state has in addition a photon remnant, keeping approximately 
the original photon direction. The two diagrams represented in Fig. are of the same 
order in the coupling constants, al^ag- In fact the resolved diagram involves the photon 
parton densities which are of order aemlf^s-, as we discussed in section |, bringing the 
power counting to the same result. 


The kinematics of hnal state particles or jets are generally specihed by their transverse 
momentum pt with respect to the beam axis (or equivalently the transverse energy Et) 
and by the pseudorapidity rj. The dehnition of a jet requires a jet algorithm at both 
the experimental and the theoretical level. We will describe in some detail the relevant 
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definitions. At the Born level the knowledge of Et {Eti = Et 2 ), Vi r ]2 for the two 
final state partons, identihed with the hnal jets, allows one to obtain the energy fractions 
x.y and Xp of the initial partons from the photon and the proton respectively. Let us write 
the four-momentum balance in the hard 2^2 subprocess: 


l/x^k-hXpP = ji -h j2 


(42) 


where k, P, ji, j 2 are respectively the four-momenta of the incoming electron and proton 
and of the two hnal partons. From this one obtains: 


E^ 


^ 2E. 


2Er 


^ (e-^1 + e-^^) 

7 

Ex 




This holds only at LO. At this level direct photon events have x^ 


x^ < 1 . 


(43) 

(44) 

1 and resolved ones 


The distinction between direct and resolved photon diagrams becomes theoretically 
ambiguous in NLO of QCD. Both components are related to each other through the 
factorization scale at the photon leg. The dependence of the NLO direct cross 
section cancels against the dependence in the resolved cross section via the photon 
structure function. In a consistent NLO calculation both the components have to be 
summed using the same scale. In general one defines observables x^^^ and Xp^^ 
which are determined from the two highest Ex jets in the event: 


OBS _ Exie Ex2e 


xr = 


2 R 


(45) 


OBS Exie'^^ + Ex2e'^^ 


x„ = 


2E^ 


(46) 


Direct and resolved photon events can be experimentally defined by a cut on x^^^ when 
two or more jets are measured in the hnal state 


The cross section for n-jet production is expressed formally as a convolution of the 
partonic cross section with the photon and proton parton densities: 


da{e e jets X) = [ dy x 

Jo 


0 

X E 

a,6 


(47) 


f dxp fi>/p{Xp,Mp) dd{a b ^ jets) 

0 


where a and b are the parton types from 7 and p, and Mp the factorization scales 
for the photon and proton parton densities. For a direct photon 0 = 7 and formally 
= S{l-x^). 
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For the production of a particle h the inclusive cross section is expressed by introducing 
the fragmentation function Dhjjyz)^ which gives the probability to obtain h from a hnal 
state parton c, with momentum fraction z = Ph/Pc- If can be written: 


X 


da{e p ^ e h X) = [ dy <F^/e(?/) x (48) 

Jo 

lo lo lo da{ab^cX) 


where Mj is the factorization scale of the fragmentation function, which is usually set to 
the p"^ of the produced particle. The fragmentation functions need to be parametrized 
from experimental results 
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The validity of equations (^7[) and (^51) rests on the QCD factorization theorems (for 
a review see [0), which state that universal parton distribution functions can be defined 
for a given hadron independently from the process. Thus comparing QCD calculations 
with experimental data may serve to constrain the parton distributions or, when the latter 
have negligible uncertainties, to test the parton dynamics. 


7.1 Inclusive charged particle distributions 


Inclusive distributions of charged particles in photoproduction have been measured by 
both HI 1^ and ZEUS |^. In Fig. ^ the transverse momentum distributions are 
shown averaged over the pseudorapidity interval —1.2 < p < 1.4 (corresponding to the 
region —3.2 < p* < —0.6 in the yp center of mass frame, in the photon hemisphere) as 
measured by ZEUS at {W) = 180 GeV for tagged photoproduction events. The events 
have been separated in a non diffractive sample and two diffractive samples with different 
average values of the diffractive system mass Mx of 5 and 10 GeV. The separation is done 
on an event by event basis, labelling as diffractive the events with a rapidity gap around 
the proton beam direction larger than a fixed cut. Soft interactions of hadrons can be 
described by thermodynamic models [Q, which predict an approximately exponential pt 
spectrum of the kind: 

1 d’^N 


Nev dpj^dp 


= exp I a 


b\JpT + ml 


(49) 


where is the pion mass. This expression is fitted to the ZEUS data in the range 
0.3 < Pt < 1.2 GeV resulting in the solid lines. The fitted values of the slope b are 
compared to results from pp and pp data as a function of the center of mass energy in 
Fig. The slope of the ZEUS non-diffractive spectrum agrees with data from hadron- 
hadron scattering at similar energies. The diffractive slopes agree better with hadronic 
data taken at a lower energy: here they are plotted at 5 and 10 GeV, the values of the 
invariant mass of the dissociated photon, in agreement with data on proton diffractive 
dissociation at hadron colliders ||^, Q . 
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The pt spectrum for non-diffractive events in Fig. ^ deviates from an exponential 
at high Pt values. This is expected from QCD as a result of hard scattering between 
constituent partons. The high pt behaviour can be approximated by a power law: 


1 <PN 
Nev dp'^dp 


A 


Pro/ 


(50) 


with pto = 0.54 GeV and n = 7.25. In Fig. ^ HERA data are compared to data from the 
WA69 m fixed target photoproduction experiment at {W) = 18 GeV. This shows that 
the Pt spectrum becomes harder as the center of mass energy increases. The figure also 
shows fits of the form (0) to pp data from UAl and GDF at several different energies 
The Pt spectrum for photoproduction at HERA is clearly harder than the one 
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from pp interactions at a similar energy, and is in fact similar to pp at s/s = 900 GeV. 


High energy yp and pp interactions show a strong similarity for the production of 
particles at low pt, while the rate of high pt particles is quite different, revealing the two¬ 
fold nature of the photon. In fact the most appealing reason for the observed difference 
is that at high pt the pointlike coupling of the photon becomes important: thus yp 
interactions are not just like vrp, instead direct and anomalous processes may occur. 
Other possible explanations could be the different pseudorapidity ranges used by the 
various experiments ( |p| < 2.5 for UAl, |p| < 1 for GDF, —3.2 < p < —0.6 for ZEUS) 
or harder parton distributions in mesons (to which the photon converts according to the 
VDM) than in baryons. 


HERA data have also been compared to a NLO QGD calculation and are in 
good agreement in the region of high pt- The high pt inclusive particle cross section 
has recently been considered an important tool to extract precise parton distributions 
in the photon in the poorly known regions of > 0.8 and < 0.05 [^. It has the 
advantage of being insensitive to possible multiple interactions or soft underlying event 
effects, discussed in section |^, as opposed to jet measurements. 


7.2 Jets 


Apart from the long tail in the pr distribution of particles, the observation of jets is a 
direct means to establish the occurence of hard scattering between constituent partons. 
A jet is a bunch of collimated particles coming from the fragmentation of a high energy 
quark or gluon. Jet measurements need a jet definition criterion or jet algorithm. Until 
now cone algorithms ||100| , |101|| have been used more widely in photoproduction at HERA, 
given the framework similar to hadronic collisions. 


Experimentally one starts from the calorimeter cells with energies above a threshold. 
For each cell the pseudorapidity p and azimuth 0 are determined from the segment joining 
the interaction vertex to the geometric center of the cell. Each cell with transverse energy 
above a few hundred MeV is considered as a possible seed for the jet Ending. The seeds 
are combined if their distance in the p — 0 space, R = v^A0^ -|- Ap^, is less than some 
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cut (usually 1 or sometimes 0.7). Then a cone of radius R is drawn around each seed and 
the energy of all the included cells is summed up in a cluster. The cluster axis is dehned 


along the Snowmass convention ||101|| from the i^T-weighted mean of the cells inside the 
cone. A new cone of radius R is then drawn around the calculated axis of the cluster 
and the axis is recalculated. The procedure is iterated until stability is reached. At the 
end it may happen that two clusters overlap. In this case if the shared transverse energy 
is more than 75% of the total Et of one of the clusters they are merged, otherwise the 
common cells are assigned to the nearest one. In the end a cluster is called a jet if its 
transverse energy exceeds a minimum value. At HERA this is set typically in the range 
of 6 - 8 GeV. 


Inclusive jet production has been measured by ZEUS ||102|| for > 8 GeV in the 


laboratory pseudorapidity range —1 < < 2, from 0.55 pb“^ of untagged data (Q^ < 

4 GeV and 0.2 < y < 0.85). Jets have been dehned with a cone radius R = 1. Similar 

from a lower statistics of tagged data 
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measurements have been reported by HI 
(290 nb-i) with < 0.01 GeV^ and 0.25 <y< 0.7. Fig. 


shows the differential cross 


section da/dE^* integrated over two rj ranges: —1 < < 2 and —1 < < 1. Data 

are corrected to the hadron level for detector effects by a complete detector simulation. 

NLO QGD calculations of inclusive jet cross sections have been carried out by several 
groups ||104|, |105|, |106| , |107| , |108|] . These calculations involve at most three partons in the 


hnal state, from which jets are dehned by the same cone algorithm as used for data. 
NLO corrections can be very large, up to a factor of two. This is rehected in a strong 
dependence of the LO results upon variations in the energy scale at which the strong 
coupling constant is evaluated or in the factorization scales and M^. In LO these 
scales are often set at about the or {pt/2)‘^ of the hard jets. For those values the size 
of the NLO corrections is about its minimum. The stability of the NLO results against 
variations of the scales is much improved. The NLO curves plotted in Fig. ^ |p.06|| are 
parton level results, not corrected for hadronization. They are consistent with the data. 


Fig. ^ presents da/dp^^^ for jets with E^^ greater than 8, 11 or 17 GeV. Gomparison 
with the theoretical NLO results ||106|| shows a preference for the GS parametrization over 
GRV in the region pEt < where the contribution from high x.y partons is important. 
This is the photon fragmentation region, where the calculation is most reliable, thus giving 
a handle to test the photon structure. In the region > 1 the data, particularly at the 
lowest E^^, are systematically above the predictions. Here the fragmentation corrections 
are sizeable and moreover the data are biased by a pedestal energy which has not been 
subtracted. This is thought to arise from a combination of soft and semi-hard interactions 
between the beam remnants in resolved photoproduction and it will be subject of a further 
discussion in the next section. A thorough understanding of this region is essential to 
carry out quantitative QGD tests in the low x.y region accessible to HERA measurements 


{Xry > 0 . 01 ). 


The largest experimental errors are due to the absolute energy scale of the calorimeters, 
whose uncertainty is within ±5% for the measurement of hadronic jets in both ZEUS and 
HI. This has been determined comparing real data and Monte Garlo simulated events. 
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HI studied the Et imbalance between the scattered electron and the hadronic jet in high 
events. ZEUS studied the ratio between (measured by the calorimeter) and the 
summed pt of the charged tracks associated with the jet (measured by the central tracking 
detector), for jets with for jets with 1 < pEt < 2 the Et imbalance has been 

studied in dijet events with the second jet in the central region. Due to the exponential 
fall of the jet cross section with increasing the 5% energy scale error translates into a 
dominant systematic uncertainty of about 20% on the normalization of the cross sections. 


This uncertainty is almost completely correlated from point to point in figures ^ and ^ 
and is indicated as a shaded band. It has to be interpreted as the maximum systematic 
shift that might be applied overall. The cross section measurements could then move all 
up or down of the indicated amount. 

It is expected from QCD that jets become narrower with increasing energy. The jet 
shape, that is the transverse energy flow versus A0 or Ap, is approximately gaussian 
around the position of the jet axis. The full width T at half maximum of a jet is expected 
to decrease as 1/Et^ ||109|| . HI has measured the dependence of T on the scaled Et of the 
jet, shown in Fig. 0' The height of the gaussian is determined after subtraction of the 
pedestal energy. The jet transverse energy is normalized to the ep center of mass energy 
|1U9|| . HI data are compared to similar measurements from jet profiles in 
The decrease of T with increasing is apparent and, within the 


according to 
pp interactions [EH 


errors, pp and pp jets show the same behaviour. 

Direct and resolved photon samples have been experimentally dehned for the hrst time 
by a cut on the variable (^) in events with at least two reconstructed jets by ZEUS 


The observed distribution after the selection cuts is shown in Fig. It is 


compared to the predictions of two Monte Carlo programs, HERWIG [GT] and PYTHIA 
im. which include LO QCD matrix elements plus parton showers and hadronization 
models. Resolved and direct photon processes, as implemented in the Monte Carlos, 
populate different x^^^ regions. The peak at high values of x^^^ is clearly attributed to 


direct photon processes. Its shift from x^^^ = 1 is due to higher order QCD corrections 
and fragmentation effects, which are simulated in the Monte Carlos. Given this clear 
separation direct and resolved photon processes can be operatively defined by a cut at 
x^^^ = 0.75. In the low x^^^ region the Monte Carlo curves fall considerably below the 
data. The reason is the energy excess observed in the forward p region in comparison 
with the Monte Carlo expectations, which is discussed in the next section. 


Dijet cross sections have been measured by ZEUS separately for x^^^ > 0.75 and 
xOBS ^ g yg from 0.55 pb“^ of untagged data |^. Jets have been required to have 
Et^ > 6 GeV. Moreover a cut on the pseudorapidity separation between the two jets has 
been applied: \Ap\ = pi — P 2 < 0.5. The cross section has been measured in terms of 
p = {pi+ P2)/‘2-, the average pseudorapidity of the jets. Assuming two jets with the same 
E^^ the kinematic variables x'?^^ and x^^^ can be expressed as: 


nSS ^ lAllcOshA 


( 51 ) 
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( 52 ) 


= 


, A?7 


For Arj = 0 the hyperbolic cosine takes its minimum value of unity and thus for given 
and fj the minimum values of and x^^^ are probed ||113|| . Moreover there is a strong 
correlation between fj and x^^^ in the direct cross section and between f] and yx^^^ in the 
resolved cross section. In Fig. da/dfj is shown for events defined as direct > 0.75), 

together with results of a LO QCD calculations using GS2 parton distributions for the 
photon and different parametrizations for the proton. The agreement is reasonable and 
improves when comparing with a QCD Monte Carlo allowing for transverse momentum of 
the incoming partons, parton showers and hadronization effects. NLO QCD calculations 
have become available only recently. The hrst results compared to ZEUS data had resolved 
diagrams calculated only at LO [p.l4| , |108|| . These calculations improved the agreement 


with data. Moreover they revealed an interesting sensitivity of the dijet cross section 
at high x/^^^ on the quark content of the photon, with data preferring CS over CRV 
parametrizations. 


The resolved cross section {x^^^ < 0.75) is shown in Fig. |^. Data are compared to a 
LO calculation using MRS (A) proton parton distributions and different parametrizations 
for the photon. Here only the shape can be reproduced by the theoretical curves, the 
normalization is below the data by a factor of 1.5 to 2 (except for LACS which is excluded 
by independent measurements from TRISTAN, as we showed at the end of section 
This discrepancy can again be attributed to underlying soft or semi-hard interactions 
between the beam remnants. 


Complete NLO calculations for dijet cross sections have been performed very recently 
|115|| and compared to preliminary ZEUS data from a larger data sample. They obtain a 
good agreement for direct photon cross sections and even for resolved photon provided a 
minimum E^r^^ >11 CeV. For jets with lower E^* the comparison is still not good and 
the possibility of multiple interactions is thus conhrmed on a more solid base. 


Besides the sensitivity to photon parton distributions, jet data can test the QCD 
dynamics of the hard parton scattering. In the dijet analysis described above, the cut 
IAt^I < 0.5 constrained the angle 6* between the jet-jet axis and the beam axis in the 
dijet center of mass system to be close to 90°. In fact one gets: 


cosO* 



( 53 ) 


Experimentally one can measure only the absolute value of cos9* because the two jets 
are indistinguishable. Different angular distributions are expected for direct and resolved 
processes. Direct processes, as the diagram in Fig. ^T|/a, have a quark propagator in the 
s, t or u channel, with t and u channel processes dominating. The dominant resolved 
processes gg —>• gg, qg —>• qg, qq qq, as the diagram in Fig. pl|/b, have t or u channel 
gluon exchange. Due to the different spin of the propagators, the angular dependence 
of the cross section is thus approximately (1 — |cos6**|)“^ for resolved processes and (1 — 
|cos6'*|)“^ for direct processes. A steeper cross section is then expected at small scattering 
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angle for resolved processes 
calculations 
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This property is expected to be preserved even in NLO 
The ZEUS experiment has measured the angular distribution in dijet 


events with Ufu > 6 GeV and \f]\ < 0.5 ||117|| . The last cut keeps events with small boost 
of the dijet system in the HERA frame, which have the jets well inside the calorimeter 
acceptance: 

V ^ Vboost = l;ln ( ) (54) 


EpXp 

X'y^'y, 


A cut on the jet-jet invariant mass is applied at Mjj > 23 GeV, to avoid the bias due 
to the cut. The resulting (icr/(i|cos6**| is shown in Fig. for direct and resolved 
events, dehned as before by the cut on . There is good agreement between data and 
theory. The different shapes are not an artefact of the cut. Monte Garlo studies 

have conhrmed that the shape is a characteristic of the spin of the propagator: even 
removing the x^^^ cut resolved photon events generated with Monte Garlo maintain the 
same shape. 


7.3 Underlying event energy 


Let us come now to the problem of the observed pedestal energy. Jet shapes are shown 
in Fig. 32/a-b for direct and resolved event samples. Here the transverse energy flow is 
plotted as a function of the pseudorapidity distance from the jet axis. In the direct case 
(^oss ^ 0.75) the jet shape is well reproduced by QGD Monte Garlos such as HERWIG 
or PYTHIA, including parton showers and hadronization. Instead in the resolved case 
^ 0.75) there is a large discrepancy on the forward side: data show a considerably 
larger amount of energy. The same effect was observed by both experiments without a 
cut on x^^^ for jets in the forward hemisphere [[I02| , |118|| . The effect can be attributed 
consistently to resolved photon events, since they are largely dominant for this conhg- 
uration. A possible explanation is the occurence of additional interactions between the 
proton and photon remnants, producing a soft underlying event or a superposition of 
soft and semi-hard multiple interactions. These phenomena should not happen in direct 
photon interactions where there is no photon remnant. 

Multiple parton scattering has been studied in pp interactions in analyses of multi-jet 
events m- Recently a strong signal for it has been reported by the GDF Gollaboration 
IM. The jet shapes and the pedestal energy in pp interactions have also successfully 


been described by a multiple interaction model ||121| . 


Multiple interaction models are currently implemented in many Monte Garlo programs 
like PYTHIA ||121|| and HERWIG ||122|| , which cover both hadron-hadron and photon- 
hadron interactions and in PHOJET ||123|| , which is specihc to photoproduction. The 
simplest model in PYTHIA generates events starting from the LO QGD jet cross section, 
where the transverse momentum cutoff is set quite low, in the range of 1 — 2 GeV. 
In this region the LO parton cross section diverges and becomes larger than the total 
non diffractive cross section and- The so-called unitarization procedure consists schemat¬ 
ically in allowing more than one interaction per event, with mean (n) = apartoniPr’''^)/^nd 
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and fluctuations calculated from Poisson statistics. The resulting distributions are quite 
sensitive to (n), the cutoff and the chosen parton densities. A similar model has 
recently been implemented in HERWIG. The PHOJET generator is intended to simulate 
at once all the components contributing to the total photoproduction cross section. It 
is based on the two component Dual Parton Model ||124|| and includes multiple soft and 
hard interactions on the basis of a unitarization scheme 
dependence on the 
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Due to this scheme the 
separating the soft and hard regions is rather weak. The soft 


part is described by soft pomeron phenomenology, with parameters tuned to pp and low 
energy yp interactions. 

These models have been used by the HI Collaboration in a detailed study of the 
properties of photoproduction events related to underlying multiple interactions m. 
The average level of Et outside jets has been studied in dijet events as a function of 
, In this study the pseudorapidity difference Ap between the two jets was required 


X. 


to be Ap < 1.2, to avoid a possible mismatch of the photon remnant as a hard jet. 
In Fig. ^ the average transverse energy density {Ej) /ApAcj) is shown as a function of 
It is determined from the summed Et in the central region of the yp collision 
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7 


a:, 

7 

(—1 < p* < 1), excluding the energy contained in cones of radius 1.3 around the jet 
axes. The long-dashed line is the energy density found in minimum bias events, where of 
course xy" is undefined. The energy density increases significantly going from the direct 
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photon {x. 
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1), where it is close to the level found in minimum bias events, to low 
events. The increase is much higher than what is expected from QCD Monte Carlos 
without multiple interactions like the PYTHIA dotted curve, and cannot be attributed to 
gluon radiation or fragmentation effects. Instead data can be described by models with 
interactions between the beam remnants as PHOJET or the PYTHIA dashed line. 


This becomes more convincing when one looks at the energy-energy correlation with 
respect to the central region. The pseudorapidity correlation function D(p*) is defined as: 


n{p*) 


1 ET,n*=0 > —ET,r)*=0,i)i< Ex^n* > —ET,r]*,i) 


(55) 


HI has measured this quantity with respect to p* = 0, where the average transverse energy 
density is at the maximum, for a sample of high Et events. The result is shown in Fig. 
data show short range correlations around p* = 0 and long range anticorrelations in the 
photon hemisphere, with a minimum at p* ~ —1.8. From PYTHIA without multiple 
interactions one would expect stronger correlation and anticorrelation than data show. 
Adding multiple interactions the correlation is reduced in both signs and agreement with 
data greatly improves. The conclusion is that the addition of uncorrelated energy to the 
events results in the correct average Et level in the underlying event (Fig. |^) and also 
gives the correct correlation strength (Fig. |3^ . 


7.4 Unfolding the gluon density in the photon 

The study of the underlying energy carried out by HI has also been used to correct their 
measured jet energies by subtracting the pedestal. The applied correction is parametrized 








as a function of and has size of 0.3 — 2.3 GeV. From an analysis of dijet events with 
pedestal subtraction HI |[126|| has unfolded LO parton variables, based on the PYTHIA 
Monte Carlo with multiple interactions included. This introduces some model dependence 
on the results. The data refer to 290 nb“^ tagged photoproduction data from the 1993 
run: the hnal sample is not large, 292 events with selected two jets with > 7 GeV and 
0 < < 2.5. First it is demonstrated the need of a gluon component in the photon and 

then the gluon momentum distribution is unfolded at LO. The distribution of the selected 
events in the unfolded is shown in Fig. |^. The full histogram represents the expected 
contribution of quarks and antiquarks in the photon, the dashed histogram the direct 
photon contribution, which is added on top of the other. These two are simply obtained 
from a LO QCD calculation at the parton level, using GRV LO parton distributions for 
the proton and the photon. The proton parton distributions at the relevant momentum 
fraction are well known and do not weaken the resolving power on the photon parton 
distributions. Data are well above the histogram at low which is right the region where 
the contribution of gluons in the photon is expected to be important. The need for a gluon 
component has been checked with a f^st on this distribution, trying to fit data with 
LO QCD predictions without gluons in the photon. Only the sum of direct-y and quark- 
antiquarks in 7 has been considered and systematic corrections have been applied to both 
data and theoretical prediction in the sense to decrease the discrepancy. The resulting 
probability is found to be 0.1%. The unfolded gluon momentum distribution is shown in 
Fig. 1^ divided by Oem- The average scale of the selected events is = 75 GeV^, 

where is used as factorization and renormalization scale for the QCD calculation. 
The resulting gluon distribution is compared to LO parametrizations of GRV-LO, LACl 
and LAC3. This latter is clearly excluded, in agreement with previous observations at 
TRISTAN m 


HERA 118 


and LEP 127 


The strong rise of LACl at low x is 
disfavoured, while GRV-LO is consistent with data. This measurement constituted the 
first extraction of a LO gluon density down to x^ = 0.04. 


7.5 Photon remnant 


As we have seen, direct and resolved processes are physically distinct as is apparent from 
the x^^^ distribution in dijet events, and enriched samples can be obtained by cutting on 
this distribution. Another possible separation criterion relies on the different topology of 
the events. Figures 37 and 38 show respectively a direct and a resolved event candidate as 
reconstructed in the ZEUS detector. In these pictures the z-axis (proton beam direction) 
points to the left. The direct candidate has no activity in the rear calorimeter, i.e. in the 
photon direction. Instead in the resolved candidate several particles are found to go in 
the photon direction. Both the events have in addition two high Ex jets. The activity in 
the rear direction is interpreted as coming from the fragmentation of the photon remnant, 
which is present only for a resolved photon. 


The ZEUS Collaboration has studied the properties of the photon remnant in events 
with two hard jets with E^^ > 6 GeV and <1.6 ||128|| . To look for a photon remnant 


jet, emerging with small px at small angles in the photon direction, the cone algorithm is 
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not well suited. The jet finding was instead carried out with the fcr-algorithm ||129|| . It is 
applied starting from single calorimeter cells, as if they were particles. The prescription 
of the algorithm is to evaluate for each pair of clusters (initially single cells) the quantity: 


k'^ = 2 min{E‘f, Ej) (1 — cosOij) 


(56) 


This is the transverse momentum of the lower energy cluster with respect to the higher 
energy cluster, for small angular separation 9ij. The proton remnant (escaping detection 
in the forward direction) is considered as a pseudoparticle with infinite momentum along 
+z. After all the values have been calculated, the two clusters with minimum k"^ are 
merged and their four-momenta added to give the four-momentum of the new cluster. The 
procedure is iterated until for all pairs the resolution variable Y becomes larger than some 
threshold Ycut- Here Y is defined as: Y = k'^/, where Eif^ is the total transverse 
energy in the event. In the ZEUS analysis this Y^ut value has been chosen to vary from 
event to event in order to find always three clusters (apart from the proton remnant), 
ordered according to their After comparison between data and Monte Carlo events 
the photon remnant has been identified in most cases as the source of the third cluster, 
which appears typically in the photon fragmentation region. The final selection required 
for it a minimum energy E^ > 2 GeV and pseudorapidity r/s < —1. After these cuts 
the distribution is found to peak at low values as expected from resolved photon 

events. In Fig. [39|/a-b the pseudorapidity and transverse momentum distributions are 
shown for the third cluster, identified as the photon remnant. A comparison with the 
standard QCD simulation by PYTHIA shows that data are shifted towards higher px and 
more central rj values. The transverse momentum for the photon remnant is generated 
in PYTHIA with gaussian distribution dN/dk^ ~ with default fco = 0.44 GeV. 

This parameter determines the hardness of the kx spectrum. A fit to the data requires 
a higher k^ of 1.7 — 1.9 GeV. This observation agrees with studies suggesting that NLO 
contributions or fluctuations of the photon into qq pairs with high virtuality may lead to a 
photon remnant with sizeable transverse momentum with respect to the photon direction 
|13CI|| . The properties of the jet from the photon remnant have been compared to those of 
the hard jets, assumed to originate from the hard parton scattering. Fig. ^/d shows the 
mean value of the particle transverse energy with respect to the cluster axis as a function 
of the cluster energy; Fig. ^/e shows the average values of the scalar sums of transverse 
and longitudinal particle energies with respect to the cluster axis as a function of the 
cluster energy. Particle quantities are measured from calorimetric islands, i.e. groups 
of adjacent calorimeter cells, and corrected to the hadron level by a complete detector 
simulation. Figure ^0|/f shows the energy flow as a function of the distance from the 
cluster axis, expressed as 1 — cos 9. In all the plots the photon remnant and the hard 
jets exhibit the same hadronization properties. The average particle transverse energy 
with respect to the jet axis {EJp) is of the order of few hundred MeV and increases only 
weakly with the jet energy. The longitudinal component of the jet energy increases much 
faster than the transverse energy w.r.t the jet axis. Fig. |40|/ a-b-c compares the photon 
remnant distributions with expectations from PYTHIA; they are in good agreement after 
the tuning of the intrinsic kx distribution of the photon remnant. 
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7.6 Open charm production 


Heavy quark photoproduction has been the subject of many theoretical studies ||131 . 


Recent NLO QCD calculations for the HERA energies have been published for both 
total and differential cross sections ||132|| . A large amount of data exist from hxed target 


experiments with center of mass energies below 20 — 30 GeV |p.33| . 

In leading order the dominant production mechanism is expected to be the direct 
photon-gluon fusion process ■yg QQ, with Q = c,b, which is represented by the diagram 
in Fig. ^/a, where the fermion line is a heavy quark. This is mostly a photoproduction 
process, with a tail at high which is not dominant. Charm production is much more 
abundant than beauty production; the predicted cross sections at HERA are about 1 /ib 
for charm and about 10 nb for beauty. Monte Carlo studies show that the dominance of 
charm production still holds after the application of experimental tagging and selection 
criteria, that is why for the present analyses, with the achieved HERA luminosity, beauty 
production has not yet been considered. In LO QCD charm production can occur also 
in resolved photon interactions, mainly through the process gg cc. The relative frac¬ 
tion of direct and resolved photon is very sensitive to the photon parton densities. For 
parametrizations not steeply rising at low x.y, like GRV, the resolved contribution is of 
the order of 20%, while for LACl it is more than 50% of the total cross section. An¬ 
other strong source of variations in the theoretical predictions is the value of the charm 
mass rric, which is generally taken in the range 1.2 — 1.8 GeV. Variation of the mass in 
this range changes the total cross section by more than a factor of 2. Charm production 
is somewhat at the boundary of applicability of perturbative QCD, due to the relative 
smallness of the charm mass, which is here the relevant scale for the perturbative expan¬ 
sion. The renormalization scale is set to rric and the factorization scales for proton and 
photon parton distributions are set to 2mc, according to ||132|| . Large variations occur by 
changing the scales. The large spread in the predictions justifies the interest in measuring 
the charm production at HERA, aiming to hx the input parameters of the calculation. 
The most interesting topic may be the determination of the gluon density in the proton, 
due to the predicted dominance of the direct photon production in the experimentally 
accessible region. In contrast to the measurements carried out from the scaling violations 
in DIS [P3| , |134|, the gluon density enters directly at the leading order of QCD in charm 


production. Several possibilities to extract it have been investigated ||135|, 


HERA experiments have measured the charm production mainly through the recon¬ 
struction of the D*(2010)^ meson decays. The golden channel is 

very small mass difference between D*^ and D^, AM = 145.42 MeV; the slow pion vr^ has 
only 40 MeV momentum in the D* rest frame. This special kinematic constraint allows 
AM to be measured more accurately than the D* mass itself | 137 . 
reconstructed in the channels K~tt^ or K~ 

decay chain from the D*~ is implicitly considered. 


IT IT TT 


The has been 
The charge conjugate 


Results on the total cross section yp —>• ccX have been published by both ZEUS 
and HI [|139||. HI results come from 1994 data, amounting to 2.77 pb“^ tagged 
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photoproduction data < 0.01 GeV^, 159 < 11^ < 242 GeV) and 1.29 pb“^ untagged 
data < 4 GeV^, 95 < hh < 268 GeV). The accepted kinematic region is pt{D*) > 
2.5 GeV and —1.5 < y{D*) < 1, where y indicates the rapidity. The pt cut depends 
on the minimum pt needed to measure the tracks from the D* decay. The rapidity cuts 
select a well understood region of the tracking detectors. No particle identihcation has 
been used yet; well measured tracks are simply combined to form invariant masses. These 
measurements are shown in Fig. 41 as a function of the average yp center of mass energy 
VF, together with lower energy measurements from fixed target experiments ||133||. The 


cross section rises by almost one order of magnitude compared to the energy range of fixed 
target experiments. The curves are NLO QGD calculations ||132|| with parton density sets 
MRS(G) ||140|| for the proton and GRV-HO for the photon. They limit the variation 
produced by a change of the renormalization scale p in the range vncl‘1 < p < 2mc- 
The limited acceptance of the HERA experiments makes the measurement of the total 
charm cross section rather model dependent, because of the need to extrapolate from 
the accepted kinematic region. The extrapolated a^p increases by 75% if LAGl photon 
parton densities are used instead of GRV-HO; it decreases by 35% if the MRS (A’) proton 
parton densities are used instead of MRS(G). These uncertainties are included in the 
error bars. Data are not able to discriminate between different available gluon density 
parametrizations, even though they favour rising distributions at low x. The visible cross 
section is instead almost insensitive to the choice of the parton distributions and to the 
mixture of direct and resolved processes assumed in the Monte Garlo to calculate the 
efficiency. Then cross sections within the experimental cuts can be safely compared with 
theory in a model independent way. 


Differential cross sections versus the transverse momentum or the (pseudo)rapidity of 
the D* have been measured by both the experiments [ 139 , |141|| . ZEUS data come from 
3 pb“^ of the 1994 run. Fig. ^ shows the differential cross section da/dpT for the observed 
The theoretical curves are two very different NLO QGD calculations: the dotted 
and the dashed-dotted curve stem from an exact fixed-order NLO calculation including 
the charm mass [|132|| ; the full curve is instead a NLO resummed calculation, where the 
charm mass is neglected in comparison to the c-quark transverse momentum [ 143 , 144 
(limit of large pT/rric)- The latter approach has been pursued also by another group 
B HHl and is sometimes referred to as the massless approach. Within the limits of 


applicability of this approach the charm is considered an active flavour in both the proton 
and the photon structure. On the contrary in the massive approach charm is generated 
from light partons and is excluded from the initial state. Both the calculations take into 
account the probability for a charm quark to fragment to D*: V{c D*~^) = 26.0 ±2.1% 
1^. The D* spectrum is obtained from the calculated c-quark spectrum convoluting it 


with the Peterson fragmentation function ||146|| , which apart from the normalization is: 


/W 


z{l -^ 

2 ; 1 — z 


-1 


( 57 ) 


Here 2 ; is the fraction of the charm quark momentum taken by the D*, and Cc = 0.06 is 


a free parameter which is htted from experimental data ||147|| . The massive calculation 
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with default parameters gives the dash-dotted curve. It has the correct shape but the 
normalization is about half that of data. With an extreme choice of rric = 1.2 GeV and 
renormalization scale equal to half the transverse mass the dotted curve is obtained. The 
effect of changing the charm mass in the range 1.2 — 1.8 GeV gives only a ±20% variation 
inside the accepted region. Instead the resummed calculation with equal parameters gives 
the correct normalization. This agreement is somewhat unexpected at least in the lowest 
Pt region covered by the data, where the massive calculation should be definitely better. 
A possible explanation to this situation is suggested in |[148|| . Here it is argued that the 
Cc parameter in the Peterson fragmentation function (^) varies sensibly depending on 
whether the fit to data is done using a LO or a NLO QGD calculation for the partonic 
piece of the cross section. Using their NLO resummed calculation the authors hnd the 
best fit for Cc — 0.015 — 0.02 to be compared to the commonly used value 0.06, which is 
found with a LO calculation. This change would increase sizeably the D* cross section 
within the experimental cuts, and would bring also the massive calculation closer to data. 
In Fig. ^ the differential cross section as a function of the D* rapidity is shown, from the 
HI analysis of untagged photoproduction data. Data are above the theoretical prediction 
of the massive calculation ||132|| in the forward region. Different curves show that the 
effect of changing the proton and photon structure functions is small. In particular using 
the LAGl instead of the GRV(HO) set for the photon makes only a small variation. 


7.7 Inelastic J/?/; production 

Another relevant process to access the gluon density in the proton is inelastic J/ijj pro¬ 
duction. Here again the leading mechanism is thought to be a photon-gluon fusion, 
''yg ^ J/%Ij g/m the so-called Golour-Singlet Model ||149|| , represented by the diagram in 
Fig. The name is due to the assumption that the cc pair giving the J/rf) lies in a 
colour-singlet state. Extraction of the gluon density from this process has been the aim of 
analyses at hxed target experiments ||150|| . There, a discrepancy was found in the absolute 
normalization, with the data being higher than the predictions by a factor from 2 to 5. 
This was ascribed to missing higher order contributions (K-factor). Recently NLO QGD 
calculations have been performed ||151| , |152|| and agreement with those data was obtained. 
HERA data have been published by HI from an integrated luminosity of 2.7 pb“^ . 
The selected decay mode was J/ip ^ The direct photon-gluon fusion process has 

to be separated by other processes with quite different characteristics, like diffractive pro¬ 
duction with proton dissociation and production via resolved photon ||153|| . They are best 
discriminated by measuring the inelasticity variable z, defined as: 


2; = 


Pyj • P 
q ■ P 


(58) 


where P, q are respectively the four-momenta of the J/ip, the initial proton and the 
photon. This expression is evaluated from measured energy and momentum variables as: 


= 


(A' Pz)lp _ (P Pz)ll> 


(P-pQ^ EhiE-Pz)h 


(59) 


33 




















with the sum extended to all the dual state particles except the scattered electron. In the 
proton rest frame z would be simply the ratio of the J/il) to the photon energy, z ~ 1 
in the elastic process yp —*• p, since in that case almost all the photon energy is 

transferred to the J/'0. ^ ~ 1 also for events where the proton diffractively dissociates. 
Experimentally if particles are detected from the proton breakup they are at very forward 
angles, with p^ ~ E, thus contributing only a little to expression ([59|). On the other 
hand resolved photon production contributes at low z, due to the reduced energy that the 
photon brings into the hard subprocess. A clean sample of photon-gluon fusion events is 
selected in the range 0.45 < z < 0.9 after detailed studies and Monte Carlo simulation 
of all the relevant processes. The cross section has been measured extrapolating to the 
low 2 ; range and is shown in Fig. |45|/a. Data are compared to a NLO calculation in the 
colour singlet model ||152|| with different gluon density parametrizations. All the curves 


assume rric = 1.4 GeV, Ejjg = 300 MeV and renormalization and factorization scales 
set to y/2mc- The agreement seems better for the distributions rising more steeply at 
low X, like MRS(G). Actually the NLO calculation is not reliable at high z and low pt 
of the J/'^. This is apparent from the differential cross section da/dp^ in Fig. ^/a, 
where the theoretical curve bends over at pf^ < 1 GeV^. A safe kinematic region to avoid 
theoretical problems is z < 0.8, p^ > 1 GeV^: this corresponds to the requirement that 
the emitted gluon is hard. The cross section in this restricted region is shown in Fig. |45|/b. 
Unfortunately the sensitivity to the input gluon distribution is decreased and all the curves 
agree with data. An exponential ht to da/dp\ with an exponential shape e~^^T gives a 
slope b = 0.39 ± 0.06 ± 0.03 GeV“^, in agreement with the theoretical calculation which 
gives b = 0.3 GeV“^ for pf^ > 1 GeV^. 


Data on inelastic J/ijj and production at the Tevatron ||154|| have shown a huge 
discrepancy with the predictions of the Golour Singlet Model. They can however be htted 
assuming additional colour octet terms in the general formalism of Bodwin, Braaten and 
Lepage ||155|| . This is not a new model: the basic idea is that the cc pair may be either in 
a colour singlet or in a colour octet state before producing the J/'0- In fact a colour octet 
state becomes a singlet simply by the emission of soft gluons. This non-perturbative 
step cannot in principle be neglected. The calculation needs phenomenological inputs 
which have to be fitted from data. To verify if HERA data are consistent with those 
observations, calculations have recently been performed to assess the effect of possible 
colour octet contributions at HERA ||156| , |157|| . The prediction would be of a very steep 


increase of the production cross section for inelastic J/^l) at high z. In Fig. ^/b da/dz 
as measured by HI is compared to the prediction of NLO colour singlet alone and to LO 
colour octet fixing its parameters to the GDF data ||154 
exclude any large colour octet contribution. 


It is clear that HERA data 


8 Hard diffraction in photoproduction 


Both ZEUS and HI have reported the measurement of a large diffractive component in 
photoproduction, about 40% of the total yp cross section |^, ^ , as discussed in section H. 
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This is expected in view of the analogy of jp to hadron-hadron interactions. Diffractive 
processes are not limited to soft interactions. The first observations of diffraction in 
hard collisions at HERA have come from DIS measurements [p.58|] . The experimental 
signature for diffractive events was the observation of a large rapidity gap around the 
proton beam direction. This is customarily specified by the variable rj^ax, defined as the 
pseudorapidity of the most forward energy deposit in the calorimeter with energy above 
400 MeV. Diffractive event samples have been defined having pmax < 1.5 — 2, which means 
that no particles are seen in the detectors within a cone of 25° — 15° degrees around the 
proton beam direction. This cut selects mainly single photon dissociation events, where 
the initial proton is not broken up. Minor components are events where the proton 
dissociates to a low mass system which is not detected or non-diffractive events with 
unusual energy flow. About 10% of the DIS events were found to have the large rapidity 
gap signature. The dependence of the event rate on W as well as the approximate scaling 
observed in the DIS data indicated a diffractive process of leading twist. The diffractive 
structure function ||159| , |16(]|| and the cross section as a function of the mass of the 

diffractive system Mx ||161|| have then been measured. Hard diffractive events have also 
been observed in photoproduction ||162| 
rapidity gap signature were found. 


Here events with high Et jets and the large 


8.1 Jets and heavy flavours from photon diffractive dissociation 


The picture of a hard diffractive event in photoproduction is shown in Fig. 47. No activity 
is seen in the forward calorimeter, although two hard jets are visible in the central region. 

The interest in hard diffractive processes is motivated by the possibility of perturbative 
QCD calculations. This can help to understand the nature of the pomeron, considering 
it as a partonic system [163|. Ingelman and Schlein |164| gave a particle-like description 


of the pomeron taking inspiration from ISR pp data ||165||. Their model was confirmed 


by the UA8 experiment ||166|| , which observed jet production in pp diffractive interactions 
with a tagged leading proton. 

Let us describe in some detail how the Ingelman-Schlein model is implemented. It 
assumes that the proton emits a pomeron: a colourless object with vacuum quantum 
numbers. The variable xjp is the fraction of the proton momentum taken by the pomeron. 
The square of the photon-pomeron center of mass energy is given by: ~ Xpys. The 

flux of pomerons from the proton is factorized in the cross section and its expression is 
derived from hadron-hadron data. It is a function oi xp and t, the momentum transfer 
from the proton. Different forms for the pomeron flux factor are currently used. The 
Donnachie-Landshoff form | 167 | is calculated in Regge theory with parameters fitted to 
hadronic data: 




(60) 


where Fi{t) is the elastic form factor of the proton, bo ~ 1.8 GeV ^ is the pomeron-quark 
coupling and a(t) = 1.08 -|- 0.25 t(GeV^) is the pomeron trajectory. The systematic un- 
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certainty related to this choice is about 30%, from the comparison of different procedures 
to extract the flux factor. 

In the Ingelman-Schlein model the pomeron is then assumed to be a source of partons 
which may interact with the photon in both direct and resolved processes, according to 
the pointlike or hadron-like behaviour of the photon. Parton densities of the pomeron 
are introduced in perfect analogy to the proton case, where P is the scaling 
variable, i.e. the fraction of the pomeron momentum carried by parton i, and fi is the 
energy scale at which the pomeron is probed. The partonic cross section is then calculated 
in LO QCD and factorization of these parton densities is assumed. If this would be the 
case they could be extracted independently of the experiment. Several forms for the 
pomeron structure have been suggested, the most extreme ones are the following: 


• iP entirely made of gluons, with a hard distribution: 

Pfg/riP,fP^) = Qp{l- P); 

• IP entirely made of gluons, with a soft distribution: 

• IP entirely made of qq pairs (only two light flavours), with a hard distribution: 

PU/AP.p') = W-P)- 


The momentum sum for the pomeron is written as: 

JO ■ 


(61) 


and in the Ingelman-Schlein model it is assumed Sjp = 1. This is however only a guess. 
Since the pomeron is not a real particle, the flux factor /jp/p and the momentum sum Sjp 
may be not independent. In the following the dependence on is neglected. 


Jet cross sections in diffractive photoproduction have been measured by ZEUS |p.6S 
from 0.55 pb“^ of untagged data. Data are corrected at the hadron level by detector 
simulation. Preselected events have a jet with > 6 GeV (at the calorimeter level) in 
the central region — 1 < < 1. This corresponds to > 8 GeV at the hadron level. 

The invariant mass of the hadronic system is reconstructed as: Mx = \JPli Ef — Pn 
where the sum is extended to all the calorimeter cells. Fig. ^/a shows the correlation 
between Mx and rjmax'- it is apparent that there is a distinct class of events with low Mx 
that have low rjmax, that is a large rapidity gap. The distribution of rjmax for events with 
Mx < 30 GeV is shown in Fig. ^/b. The simulation of non-diffractive processes by the 
PYTHIA Monte Garlo is not able to describe the tail at low rjmax- Instead the POMPYT 


m 


Monte Garlo ||169|| , based on the Ingelman-Schlein model ||164|| , describes well the shape 
this region. Fig. ^/c shows the Mx distribution after the cut rjmax < 1-8 and Fig. |4^/d 
the W distribution for the selected events. The flat W behaviour is in agreement with 
the expectations of POMPYT. 
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Figure ^ shows the differential jet cross section versus for the large rapidity 
gap events. Data include a contribution from double dissociation, where the proton 
diffractively dissociate to a hadronic system of mass M^v < 4 GeV which is not detected. 
The dominant errors are statistical, while the largest systematic error 20%) is due 
to the uncertainty on the absolute energy scale of the calorimeter. LO QCD predictions 
obtained by Monte Carlo are superimposed. Large rapidity gaps may occur in normal 
non-diffractive events as fluctuations in the pseudorapidity distribution of the final state 
hadrons. They are however exponentially suppressed with increasing gap width. The 
size of the non-diffractive contribution has been estimated with PYTHIA using MRS(D-) 
parton distributions for the proton ||184 | and GRV(HO) for the photon. The result is 
considerably lower than data and differs in shape. Data are instead in good agreement 
with the expectations of POMPYT assuming a hard gluon density inside the pomeron, 
while a pomeron made of hard quarks gives the correct shape but too low a normalization. 
A soft gluon density gives inconsistent results on both shape and normalization. These 
data are however insensitive to an additional soft component in the pomeron, present 
together with the hard component. 

The diffractive jet data are sensitive to both the quark and the gluon content in the 
pomeron. The recent HERA measurements of the diffractive structure function in DIS 
have provided information on the quark content. This appears to have both 


159, 160 


a hard and a soft component. Moreover these measurements do not favour a pomeron 
made only of quarks and antiquarks. If factorization of the pomeron parton densities 
holds, the pomeron structure would appear the same whether it is probed in DIS or in 
photoproduction. Therefore, assuming factorization, the gluon content of the pomeron 
can be studied combining DIS and ■yp jet data. For that purpose a pomeron composed of 
both hard gluon and quark distributions has been assumed, naming Cg and Cq = 1 — Cg the 
relative gluon and quark fractions. Gontributions from a possible soft parton component 
have been neglected, since the jet data are sensitive only to /? > 0.3. The distribution of 
da/dr]^^* was corrected subtracting the non-diffractive and the double dissociation con¬ 
tributions. The hrst one was estimated with PYTHIA, which gives a good description 
of non-diffractive jet production in the central rapidity region. The latter was estimated 
from data to be 15 ±10% and it has been assumed independent from according to the 
theoretical expectations. Then a fif to the measured da/dp^'^^ was carried out varying 
the gluon fraction Cg and leaving free the momentum sum (1^ in POMPYT. The results 
are shown in Fig. 0 by the solid line. The shaded band gives the Icr contour. 

From the measurement of Q^,X]p) in DIS we can find the quark contribution 

to the momentum sum, J2iPq, integrating over (3 and Xp: 




/Xip„ 


r*l 


dxp I dl3 Q^, xp) — kf ■ Y^pq ■ Ifiux 

Jo 


(62) 


where Ifiux is the integral of the pomeron flux factor over t and over the same region 
in xp and kf is a number depending on the assumed number of flavours (5/18 for two 
flavours, 2/9 for three flavours). The integral on the left hand side is performed using 
ZEUS data in the range 6.3 ■ 10“'^ < xp < 10“^. The flux factor Ifiux is calculated 
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using expression (|60D. The result is: Sjpg = 0.32 ± 0.05 (0.40 ± 0.07) for two (three) 
flavours. The diffractive structure function has been parametrized from ZEUS data 

in the (3 region 0.1 < /? < 0.8. Here the same form is assumed to hold in the whole 
range 0 < /? < 1. is approximately independent of in the range of the DIS 

measurements: 8 GeV^ < < 100 GeV^. Hence the dependence of Sjpg on is 

neglected. Moreover the values are comparable to the values of (Uf®*)^ of the selected 
jet events. Thus the DIS results are expressed as a constraint in the Ejp — plane of 
hgure ^ (dot-dashed lines): 


Ejp ■ (1 — Cg) = Ejpq = 0.32 (0.40) for 2 (3) flavours 


(63) 


Gombining the two results and considering all the systematic uncertainties one obtains: 
0.4 < Ejp < 1.6 and 0.3 < Cg < 0.8. From this nothing can be safely concluded on the 
validity of a momentum sum rule for the pomeron. However, assuming the same pomeron 
flux factor in DIS and photoproduction, this determination of the gluon content of the 
pomeron is independent on the normalization and even on the validity of a momentum 
sum rule. From this analysis one concludes that between 30% and 80% of the momentum 
of the pomeron carried by partons is due to hard gluons. 


Further information on the partonic structure of the pomeron can be obtained from 
the diffractive production of heavy flavours. The diffractive production oi J/f) has been 
discussed in section ^ because of the similarity with the light vector meson analyses. It 
illuminates the realm of hard diffraction, exhibiting a strong difference from light vector 
mesons, particularly in the steep dependence of the cross section on the center of mass 
energy. Glean data samples |^, with good statistics have already been studied at 
HERA and succesfully compared to calculations based on perturbative QGD |8^ . However 
for J/fj elastic production the theory retains some weakness related to long distance 
contributions in the cross section. Instead the diffractive production of open charm is 
in principle free from these uncertainties and recently has attracted much theoretical 

|7g 


interest m- This process has been searched for by HI ||139|| . They observed evidence 
of diffractive charm production by inspecting the rj^ax distribution, shown in Fig. ^ for 
selected D* events. Here statistical subtraction of the background under the mass peak has 
been applied. The QGD Monte Garlo representing non-diffractive production is unable to 
reproduce the distribution, while the hard diffraction model implemented in the RAPGAP 
generator ||171|| reproduces well the shape when added to PYTHIA. HI sets a lower limit 
on the diffractive cross section for pt{D*) > 2.5 GeV, —1.5 < y{D*) < 1 and rjmax < 2 in 
untagged events (Q^ < 4 GeV^ and 0.1 < ?/ < 0.8) to: a{ep —> D*^X) > 145 pb at 90% 
conhdence level. Gomparing this limit to predictions of RAPGAP, a quark dominated 
pomeron is disfavoured. In fact, the predicted cross section for a pomeron made of hard 
gluons is ~ 800 pb, while for a quark dominated pomeron it is only 30 pb. 


8.2 Hard colour singlet exchange 

Another kind of hard diffractive process leads to the appearance of large rapidity gaps 
between two high Ft jets. These events would result from the exchange of a colour singlet 
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propagator in the t-channel, as in the diagram of Fig. [^/a. The propagator conld be an 
electroweak gange boson or a strongly interacting object. The normal QCD jet production 
is instead described as a strong process with the exchange of a quark or a gluon propagator, 
as in Fig. ^/b. The colour charge of the propagator generates the colour connections 
between the hard jets. For this reason the pattern of soft gluon radiation is distributed 
over all the phase space, particularly in the central rapidity region between the two jets. 
Instead, for colour singlet exchange, the hard jets would be colour connected only to the 
beam remnants. In this case the central region would be depleted from particle production 
and could present the rapidity gap signature. This was first suggested by Bjorken ini. 


who also calculated the expected rate for this kind of events in pp interactions at the 
Tevatron {^/s = 1.8 TeV). He estimated the contribution from the exchange of two gluons 
in a colour singlet state to be 10% of the cross section for single gluon exchange. Other 
calculations valid for high t have been published based on perturbative QCD ||173| , |174 


The topology of the relevant events is shown in Fig. ^/c. Two jets are shown in the 
rf-cj) plane, back-to-back in azimuth and separated by an interval Arj in pseudorapidity. 
This is the pseudorapidity difference between the axes of the two jets. No particles are 
found between the jets in rapidity gap events. 

The gap-fraction f{Ar]) is defined as the ratio of the number of dijet events having 
a gap Arj to the total number of dijet events with jet-jet separation equal to Arj. It 
is expected to fall exponentially with increasing Arj for non-diffractive processes. This 
follows assuming a constant pseudorapidity distribution of particles between the jets and 

13 


and Monte Carlo results. Instead for 


is supported by analytic QCD calculations 
colour singlet exchange the gap fraction is expected to be almost flat with Arj \\172[ |174 


At large enough Arj the colour-singlet exchange would then be dominating, as is shown 
in Fig. ^/d. 

Allowing secondary interactions between the photon and proton remnants, the gap 
could be sometimes filled. A survival probability V has been defined in ||172|| , such that 
the observed gap-fraction is the product of V and the gap probability at the parton level, 
which is a calculable quantity. Estimates of the survival probability for pp collisions at 
the Tevatron are in the range 5 — 30% ||172|, |175||. The same quantity at HERA could be 


higher due to the substitution of the p with the 7 and the lower center of mass energy. 


Both the Tevatron experiments DO | 176|| and CDF [1^] reported the observation of 
an excess of events with two hard jets separated by a rapidity gap, with respect to the 
expectations from normal QCD processes with colour exchange. The measured excess is 
about 1% for DO, 0.8% for CDF. 


ZEUS published results based on 2.6 pb“^ of 1994 data ||178|| , requiring two jets with 
> 6 GeV, 77 -^®* < 2.5 and small boost \p\ = {rji -|-? 72|/2 < 0.75. A cone algorithm with 
R = 1 was used and the distance between the two jets was required to be | Ar/I > 2 to have 
non-overlapping cones in p. Gap events were found by inspecting the calorimeter pattern 
of energy, identifying particles by “islands” of neighbouring calorimeter cells with summed 
Et above a threshold of 250 MeV. A gap event was dehned as an event having no particles 
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between the two jets with Et > 300 MeV (after correction of the calorimeter energies). 
The result is shown in Fig. where the excess of gap events is evident in the last bin. 
The largest systematic uncertainty in data comes from the cut on the of the island and 
is included in the error bars. Data have been htted with a function /{Arj) = + /?, 

with a and (3 free parameters and k a normalization hxed by the obvious condition 
/{Ar]) = 1 at Ar] = 2. The result is shown as the solid curve, and = 1-2 for two 
degrees of freedom. The htted parameters are: a = —2.7 ± 0.3{stat.) ± 0.1(s?/st.), (3 = 
0.07±0.02(stat.)lo;o2(-s?/st.). The size of the excess has been also checked comparing data 
with the expectations from PYTHIA, simulating non-diffractive jet production. From 
this comparison the excess is found to be 7%, in agreement with the result of the ht. The 
actual number of events due to colour singlet exchange could be even higher than 7%, 
accounting for a survival probability for the gaps. The larger excess of gap events found 
at HERA in comparison to the Tevatron results could come from the same underlying 
dynamics, given different survival probabilities for gaps in the two environments. At 
HERA the contribution to the gap events from or exchange cannot account for 
the observed excess, since < 7 ■ 10“^ in this kinematic range. Instead data 

can be explained by a hard diffractive process based on the exchange of a colour singlet 
combination of gluons. 


9 Outlook 


In this paper we have reviewed the published HERA results on photoproduction. The 
understanding of the hadronic structure of the photon as well as the understanding of 
photon induced reactions has greatly improved in the last few years. HERA data have 
played a central role. This progress will continue in the next years. We will not summarize 
here the covered contents, but rather prefer to have a glimpse on the future prospects for 
selected topics. 

Most of the HERA measurements will have signihcant improvements in the next years 
with the expected luminosity increase. The hnal HERA results presented here include 
data up to the 1994 run, with integrated luminosity of about 3 pb“^. The achieved 
luminosity during 1995-96 runs is about 20 pb“^ and almost the same should be achieved 
during 1997. Future upgrade plans aim at an integrated luminosity exceeding 100 pb“^ 
per year usable by the experiments. 

Hard photoproduction processes at HERA are complementary to the ey deep inelastic 
scattering, being directly sensitive to both the quark and the gluon component of the 
photon. Data on jet photoproduction are presently consistent with the expectations of 
perturbative QCD. Ongoing analyses are studying extended domains in Ei^^ and with 
different jet algorithms. The most severe experimental limitation on jet measurements is 
the energy scale uncertainty of the calorimeters, which can however be improved to about 
2 — 3% from the 5% quoted in the results shown here. With such improvement large 
statistic samples with high E-r jets will allow more precise QCD tests. The measurement of 


40 



high pt inclusive charged particle spectra has been recently proposed as a valid alternative 
to have access to the parton distributions in the photon |^. It has the advantage of being 
insensitive to the energy scale of the calorimeter, at the price of an uncertainty due to 
the fragmentation functions, needed to describe the transition from partons to hadrons. 
The estimated accuracy on the extraction of a NLO gluon density by this method is 
about 10%. In parallel to HERA, the knowledge of the parton distributions in the photon 
(particularly the quark component) will be improved by high luminosity data from LEP2. 
From it high precision measurements of expected up to of few hundred GeV^ 

and minimum about 10“^. 


The heavy flavour studies at HERA are presently all limited by statistics. The compar¬ 
ison of open charm {D*) data with the existing NLO QCD calculations has proved to be 
an interesting and intriguing benchmark. In particular the discrepancy of the exact mas¬ 
sive NLO calculation with HERA data has to be clarified. It could be due to uncertainties 
of the charm fragmentation function or, if nothing else, it might reveal the importance of 
higher order contributions (NNLO), which are far from being calculated. Apart from the 
normalization of the cross section, it will be interesting to see if the reported excess of 
D* data in the forward pseudorapidity region will be confirmed. At reasonably high pt 
values a safe comparison of the NLO massive calculation with the calculation resumming 
the large logarithms of (the “massless” approach) should also help to understand 

the situation. Open charm data could also be used for a direct extraction of a NLO gluon 
density in the proton when a large data sample (~ 100 pb“^) will be available, as it has 

ra 


been proposed in 


Charm tagging capabilities will benefit from detector upgrades 


of the interaction region with the introduction of a silicon microvertex detector. This has 
already been installed in HI and is currently being developed in ZEUS. 


Elastic vector meson production at HERA is being studied with larger statistics while 
varying the involved energy scales t). The early results suggested the picture of a 
hard production mechanism, when a hard scale is present in the interaction, be it a high 
mass (as the J/'0 mass) or high Q'^ or t. Recent fixed target data on p production at high 
make difficult to draw safely this conclusion for light vector mesons. Detailed studies 
varying the scales will possibly shed some light on the transition from a non-perturbative 
description of the pomeron (the soft pomeron), to the hard perturbative pomeron. Hard 
diffraction data from HERA (in which the final state has jets or heavy flavour particles) 
are presently indicating that the gluon content of the diffractive exchange is substantial. 
Many different models are appearing on the market and are being tested, trying to unveil 
the nature of the hard pomeron. 


In conclusion, the study of yp interactions at HERA has yielded a lot of results. Several 
detector upgrades have already been performed (or are being done), which will result in 
improved purities of the selected samples as well as a wider kinematic coverage. Adding 
to this the expected much larger statistics, there is a bright future for photoproduction 
physics at HERA. 
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a{-fp 

VY), PD 

8.7 

± 

1.5 

± 

1.5 

± 

3.0 

3.6 

cr(7p 

Vp), EL 

17.1 

± 

1.6 

± 

3.7 

± 

1.4 

4.3 

EL + GD + PD + DD 

69.2 

± 

3.4 

± 

8.8 

± 

9.3 

13.2 


ND 

96.1 

± 

3.5 

± 

14.7 

± 

9.6 

17.9 

Total 

165.3 

± 

2.3 

± 

10.9 

± 

1.3 

11.2 


Table 1: Total and partial 'jp cross sections as measured by HI |^|. The contribution of the 
double dissociation is assumed in the range: 0 < cjdd < 40 pb. The first error is statistical, the 
second one is systematic, and the third error is the systematic uncertainty due to the assumption 
on a^D- The full error is obtained as the quadratic sum of all of these. 


Reaction 

Gross Sections (/ib) 

Data 

CKMT 

SaS 

GLM 

a{'yp Vp), EL 

17 ± 4 

17 

16 

17 

a(7p ^ Xp), GD 

26 ± 5 

25 

13 

18 

a(7p^ VY), PD 

9± 2 

7 

10 

15 

a{-fp XY), DD 

15 

15 

13 

15 


Table 2: Comparison of diffractive cross section calculations [ 
ments |Q. The data are given for a fixed value ctdd = 15/rb. 


^ 61] with the HI measure- 


reaction 

reference 

{W) [GeV] 

b [GeV-^j 

t -range [GeV^j 

7P —>• p^p 


m 

70 

10.4 ± 0.6 ± 1.1 

< 0.15 

IP P^P 


11 

70 

9.8 ± 0.8 ± 1.1 

0.073 - 0.4 

IP P^P 


M 

55 

10.9 ± 2.4 ± 1.1 

0.025 - 0.25 

7P —>• ujp 


p9| 

80 

10.0 ± 1.2 ± 1.3 

< 0.6 

-yp —>• (fp 


|70| 

70 

7.3 ± 1.0 ± 0.8 

0.1 -0.5 


Table 3: Summary of the determinations of the b slope for elastic light vector meson photopro¬ 
duction at HERA. The values of b have been obtained from a fit of the differential cross section 
of the form = A ■ in the given t range. 
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Figure 1: The kinematics of a single-tag inclusive 77 event. 



Figure 2: Quark-Parton Model (box) contribution to the photon structure function. 
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Figure 3: Recent data on F 2 {x,Q^) from TRISTAN |]^ and LEP |^, compared 


to the LO fits to all previous data of LAC and GRV Compilation from ||179| 
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Figure 4: evolution of Fj 


'2 averaged over 0.3 < x < 0.8 (from [ 180|| ). The lines are 
predictions of the FKP model with different phenomenological values of the cutoff 
parameter 0.1 (dotted), 0.5 (solid) and 1.0 GeV (dashed). The VDM contribution is 
represented by the dot-dashed line. 
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Figure 5: Parametrizations of the u-quark distribution in the photon at LO ^ 

and NLO 


3l|, from ||179|]. The NLO results are given in the MS scheme. 




Figure 6: Parametrizations of the gluon distribution in the photon at LO 
and NLO 


27, 31, 29, 37 


3l|, H, H, from 1^; 
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Figure 7: Typical leading order diagrams contributing to jet production in two-photon 
exchange processes at colliders: (a) direct photon; (b) single-resolved photon; (c) 

double-resolved photon 
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Pt (GeV) 


ee--> ee + jet + jet + X TOPAZ 



Figure 8: Inclusive (a) jet and (b) two-jet cross sections as a function of the jet transverse 
momentum for the central pseudorapidity region < 0.7 from [^. The contribution of 
the direct process only is represented by the dashed line. The other curves are predictions 
including direct and resolved processes, from different LO parametrizations ® Ezrni- 
The expectation from the DG parametrization excluding the contribution of gluons from 
the photon is also shown. 



















Figure 10; Schematic view of the central part of the ZEUS detector. 
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T1 Beam pipe and beam magnets 
"J] Central tracking de\’ice 
^ Forward tracking dje\'ice 
^ Electromagnetic LAr calorimeter 
Hadronic LAr calorimeter 
^ Superconducting coil (1.16 T) 

Y] Compensating magnet 
^ Helium supply for 


Muon chambers 
1101 Instrumented iron jxtke 
1111 Forward muon toroid 
1121 Backw. electromagn. calorimeter (BEMC) 
1131 riTJ C calorimeter 
1141 Concrete shielding 
1151 liquid argon crj’ostat 


Figure 11: Schematic view of the central part of the HI detector. 
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Figure 12: Integrated luminosity delivered by HERA as a function of the date of the year 
in the period 1992-1995. 
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Figure 13: Representative diagrams for the different classes of photoproduction processes: 
(a) elastic vector meson production 'yp ^ Vp] (b) diffractive photon dissociation 7 p — 
Xp] (c) diffractive vector meson production with proton dissociation 7 p — VY ; (d) double 
diffractive dissociation 'jp XY ; (e) non-diffractive inelastic scattering jp H. 
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Figure 14: Total photoproduction cross section measurements as function of the 7 p 
centre of mass energy Data are shown together with parametrizations which did not 
include HERA measurements: the solid curve is DL ||^ with default pomeron intercept 
ttip(O) = 1.08, while the dotted one is DL with Q:jp(0) = 1.11 (which is obtained including 
in the fit the CDF result @]); the dashed curve is ALLM 
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Figure 15: Differential cross section da/d\t\ for elastic photoproduction, 7 p — 
at {W) = 73 GeV, as obtained by ZEUS through direct measurement of the scattered 
proton. The line is the exponential £t described in the text. 



Figure 16: Exponential slope parameter 5 as a function of the photon energy in the 
rest frame of the proton. The curve is a prediction based on pomeron exchange . 
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Figure 17: Elastic vector meson photoproduction cross sections for a;, 0, J/V' as 
function of the 7 p center of mass energy. For comparison the total photoproduction 
cross section measurements are shown in the upper part of the plot. The data points at 
IF > 40 GeV are HERA measurements, the lower energy points come from hxed target 
experiments. The dashed curves are expectations from the soft pomeron model 
arbitrarily normalized at hxed target points in the elastic cross sections. 
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ZEUS 1994 




Figure 18: The elastic J/'ll) photoproduction cross section as a function of the 7 p center of 
mass energy W. In (a) the ZEUS data alone are htted to a dependence with 5 = 0.92; 
the dashed line shows the prediction of a soft pomeron model in which 5 ~ 0.22. 
In (b) HERA data are shown together with results from hxed target experiments. The 
full line is the result of a perturbative QCD calculation [Q using MRSA’ proton parton 
distributions, the dotted and the dash-dotted lines are results of two pomeron models 

IT! 
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Figure 19: Diagram representing elastic J/ijj production according to QCD-inspired mod¬ 
els with the exchange of a gluon ladder. 



|t| (GeV") 


Figure 20: The differential cross section da/d\t\ for the elastic process 7 p ^ J/ip p 
as measured by ZEUS in the kinematic range 40 < IF < 140 GeV. The result of the 
exponential fit in the range |t| < 1 GeV^ is shown as the solid line. 
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Figure 22: Inclusive transverse momentum distributions of charged particles in pho¬ 
toproduction events at {W) = 180 GeV averaged over the pseudorapidity interval 
— 1.2 < rj < 1.4. Solid lines indicate hts of equation (^51) to the data in the region 
of pt < 1.2 GeV. The dotted line shows a power law formula (^) htted to the non- 
diffractive data for p-p > 1.2 GeV. The diffractive data are scaled down by two orders of 
magnitude. 
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Figure 23: Inverse slope of the exponential fit (|4^ as a function of the center of mass 
energy ^/s for ZEUS jp data and hadron-hadron data from fixed target |p.82|| , ISR ||183| 
and UAl The ZEUS non-diffractive point is given at the jp c.m.s. energy, while 


the diffractive points are plotted at the energies corresponding to the mean value of the 
invariant mass (Mx) of the dissociated photon system. The dashed line is a parabola in 
log{s) and was fitted to all the hadron-hadron points to indicate the trend of the data. 
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Figure 24: Transverse momentum spectra of charged particles in different experiments: 
ZEUS and HI photoproduction data are compared to data from OMEGA |^, UAl 
and CDF 
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Et'" [GeV] 


Figure 25: Differential ep cross section da/dE^^ for inclusive jet production integrated 
over two ranges: — 1 < < 2 and — 1 < < 1 (scaled down by an order of 

magnitude). The kinematic range is < 4 GeV^ and 0.2 < y < 0.85. The shaded 
bands give the uncertainty due to the energy scale of the jets. The superimposed lines are 
results of a NLO QCD calculation ||106|| using MRS(D-) ||184|| proton parton distributions 
and two different photon parton distributions: GRV-HO (solid line) and GS-HO [pS 
(dashed line). 
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Figure 26: Differential ep cross section da/dp^^^ for inclusive jet production integrated 
over from three different thresholds {Et > 8, 11 and 17 GeV) in the kinematic range 
< 4 GeV^ and 0.2 < y < 0.85. The shaded bands give the uncertainty due to the 
energy scale of the jets. The superimposed lines are NLO QGD results as in the previous 
hgure. 
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Figure 27: Jet width measurements as function of the scaled jet transverse energy 
/y/s from HI ||103|| and UAl |[I10|| . HI data are normalized to the ep center of 
mass energy. F is dehned as the full width at half maximum above the pedestal energy 
determined from the Et azimuthal prohle distribution around the jet axis in HI data and 
from the pseudorapidity prohle in UAl data. The line represents a 1/El^^ £t to the HI 
data. 
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Figure 28: The uncorrected distribution (indicated with from events with 

at least two reconstructed jets measured by ZEUS. The solid and dashed histograms 
represent respectively the expectations from PYTHIA and HERWIG Monte Carlo and 
have been normalized to fit the direct peak in the data. The shaded histogram shows the 
direct contribution to the HERWIG distribution. 
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Figure 29: Differential ep cross section da/dp for inclusive dijet production in the direct 
photon sample dehned by the cut > 0.75, as measured by ZEUS. Jets are required 
to have > 6 GeV and IAt^I < 0.5 and the events to have < 4 GeV^ and 0.2 < 
y < 0.8. The shaded band gives the uncertainty due to the energy scale of the jets. The 
data are compared to analytic LO QGD calculations ||113 | using the GS2 photon parton 
distributions and different sets for the proton: GRV-LO [3^, GTEQ2M ||186|| 
m and MRSDO’ 


MRSA 
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Figure 30: Differential ep cross section da/dp for inclusive dijet production in the resolved 
photon sample defined by the cut < 0.75, as measured by ZEUS. Jets are required to 
have > 6 GeV and \ Ap\ < 0.5 and the events to have < 4 GeV^ and 0.2 < y < 0.8. 
The shaded band gives the uncertainty due to the energy scale of the jets. The curves are 
LO QGD calculations using the MRSA | 140|| parton distribution set for the proton and 
different sets for the photon: LAG1-LAG3 2^, GS2 |^, GRV [^| and DG . 
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|COS0*| 


Figure 31: Differential ep cross sections (icr/(i|cos6'*| for dijet production in direct and 
resolved pliotoproduction as measured by ZEUS, normalized to unity at Icos^*! = 0. 
Events are selected with Mjj > 23 GeV, > 6 GeV, p < 0.5,< 4 GeV^ and 


0.2 < y < 0.8. Data are compared to LO and partial NLO QGD calculations [p.l6| 
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Figure 32: Uncorrected transverse energy flow 1/N dEx/ddr] around the jet axis, for 
cells within one radian in 0 of the jet axis for (a) direct and (b) resolved photon events, 
from the ZEUS analysis in The solid and dashed lines represent respectively the 


distribution from PYTHIA and HERWIG. 
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Figure 33: Average transverse energy density (per unit area in the rj — (j) space) in the 
central region |? 7 *| < 1 outside the two jets with the highest as a function of x.y. The 
inner error bars are statistical, the outer ones result from the quadratic sum of statistical 
and systematic errors. The long-dashed line indicates the energy density measured in 
minimum bias events. The histograms show the expectations of different QCD generators 
with interactions of the beam remnants (full=PHOJET, short-dashed=PYTHIA), and 
without them (dotted=PYTHIA). 


81 















o 


* 

S3 

a 


0 


-1 


-3 


^4 


data 

PYTHIA mia 
PYTHIA 


-2 


-1 


0 






Figure 34: Energy-energy correlation with respect to r;* = 0 as a function of rj* for a high 
Et event sample ||103|| . The dashed (dotted) histogram represents the prediction of the 


QCD generator PYTHIA with (without) interactions of the beam remnants. 
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Figure 35: Distribution of the unfolded from a two-jet event sample ||126|| . Only 


the statistical errors are shown. The histograms are PYTHIA Monte Carlo calculations 
showing the contribution of quarks from the resolved photon (solid) and of direct photon 
events (dashed). 
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Figure 36: The gluon momentum distribution of the photon divided by the hne structure 
constant a as obtained in the LO analysis of HI 
compared to LO parametrizations: GRV (full) [^] . 

m- 


|126|| at a scale {pi^ = 75 GeV^, 


LACl (dashed) and LAC3 (dotted) 
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Figure 37: Direct photoproduction of jets as seen in the ZEUS detector: the figure on 
the left is a longitudinal section of the detector, in which the electron beam comes from 
the left and the proton beam from the right. Two central jets are visible, balancing each 
other in px- The upper right picture is a transverse view of the same event and the lower 
right hgure is a lego-plot showing the transverse energy deposited in the calorimeter in 
bins of f] — (p. 


84 










































ZEUS 1993 


T> 

'Z. 

TJ 0,3 


% 


0.2 - 


0.1 


0 


- 

• Data 

a 


PYTHIA 


— 

PYTHIA High k, 


- 

jj-f:. 

r-n 


- 


_L_ 

_ 1 _1_ 1 _ 1 _ 1 _1 _L_ 

_1_1_ 


-4 


-2 


^ 0.6 

V 

Q> 

P 0.4 

Q. 

T5 


TJ 

5 


0 




^3 




Data 

PYTHIA 

PYTHIA High k, 




i 



4 6 

P«(GeV) 


Figure 39: Corrected distributions of the third cluster in photoproduction events with 
two other jets with > 6 GeV and < 1.6 as measured by ZEUS ||128||: (a) pseu¬ 


dorapidity; (b) transverse momentum. The histograms are expectations from PYTHIA 
Monte Carlo: the full line is the default version, the dotted line is obtained with a higher 
intrinsic transverse momentum (kt) for the partons in the photon, as explained in the 
text. 
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Figure 40: Fragmentation properties of the photon remnant compared to those of the hard 
jets in selected resolved photon events measured by ZEUS ||128|| and to the expectations 
of PYTHIA Monte Carlo: (a) and (d) show the mean value of the average energy 

transverse to the cluster axis per particle, as a function of the cluster energy; (b) and 
(e) show the average values of the total transverse and total longitudinal (EiEl) 

energy; (c) and (f) show the energy flow around the cluster axis. 
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Figure 41: Total charm photoproduction cross section as a function of W^p-. measurements 
from ZEUS, HI and previous hxed target experiments are shown. The error bars give 
the statistical and systematic errors added in quadrature, including the uncertainty in 
the extrapolation. The lines are predictions of a NLO QCD calculation |[132|| showing the 
effect of changing the renormalization scale from /i = 0.5 (upper curve) to fi = 2 rric 
(lower curve). The parton densities used in the calculation are MRSG for the proton and 
GRV-HO for the photon. 

























Figure 42: The differential ep cross section da/dp!^ for D* photoproduction in the kine¬ 
matic range < A GeV^ and 115 <W< 280 GeV integrated over the D* pseudorapidity 
range —1.5 < < 1.0. NLO QGD predictions are shown: the dot-dashed curve is the 

result of the massive calculation of [p.32|| using MRSG and GRV-HO parton densities for 
the proton and photon respectively, rric = 1.5 GeV, renormalization scale pr = rriT and 
fragmentation parameter Cc = 0.06; the dotted curve is the same but with pR = 0.5 tut 
and rUc = 1.2 GeV; the full curve comes from the resummed calculation in ||144|| with the 
same parameters as the dot-dashed curve except the proton parton densities GTEQ4M 
I87i . 
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Figure 43: Differential 7 p cross section for D* photoproduction with respect to the D* 
rapidity y integrated over pt{D*) > 2.5 GeV. The histograms are predictions of a NLO 
QCD calculation ||132|| with rric = 1.5 GeV and different parton densities: MRSG + GRV- 
HO (solid), MRSA’ + GRV-HO (dashed), MRSG + LAGl (dotted) for the proton and 
the photon respectively. 



Figure 44: Photon-gluon fusion diagram for inelastic J/ijj production. 
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Figure 45: Total cross section for inelastic J/ijj ptiotoproduction as a function of the 7 p 
center of mass energy with: (a) 2 ; < 0.9; (b) 2 ; < 0.8 and > 1 GeV^. The curves are 
NLO QCD calculations ||152|| with different input gluon distributions and contain a 15% 


correction accounting for -0' background. 
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Figure 46: Differential 7 p cross sections for inelastic J/ip production as measured by HI 
in the energy range 30 < hF < 150 GeV: (a) with respect to 'Pr^J/'d)) integrated over 
2 ; < 0.9; (b) with respect to z, integrated over > 1 GeV^. The dashed curves are 
results of a NLO QGD calculation ||152|| in the Golour-Singlet Model, with GRV parton 
densities. The full line in (a) is a fit to the data at pf. < 5 GeV^. The solid curve in (b) 


is the expected contribution of colour octet terms in a LO calculation ||156 
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Figure 47: Dijet photoproduction with a large rapidity gap around the proton direction 
as seen in the HI detector. The left picture is a longitudinal section in which the proton 
beam comes from the right, the upper right picture is a transverse view and the lower 
right one is a lego-plot of the calorimeter energy pattern. This event has the scattered 
electron detected in the low-angle tagger, represented in the upper left corner. 
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ZEUS 1993 
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Figure 48: (a) The scatter plot of versus 77 “^^ for a jet event sample with > 

6 GeV and — 1 < <1; (b) the distribution of for the events with < 30 GeV 

along with the predictions of PYTHIA (shaded area) and POMPYT with a hard gluon 
density in the pomeron (solid line). The predictions are normalised to the number of data 
events above and below = 2.5, respectively; (c) the distribution in for the events 
with <1.8 together with the prediction of POMPYT with a hard gluon density in 
the pomeron (solid line) normalised to the number of data events; (d) the distribution in 
ppcaz events with <1.8 and the prediction of POMPYT as in (c). 


94 






ZEUS 1993 


JD 

CL 



U 


10 ^ 


10 


1 


® ZEUS Data 




soft gluon 


-1 -0.5 0 0.5 1 



Figure 49: Differential ep cross section da/< 1.8) for inclusive jet production 
with > 8 GeV in the kinematic region < A GeV^ and 0.2 < y < 0.85. The 
measurements are not corrected for the contributions from non-diffractive processes and 
double dissociation. The shaded band displays the uncertainty due to the energy scale of 
the jets. The full curves are predictions of POMPYT for single diffractive jet production 
using the DL flux factor with different parametrizations of the pomeron parton densities 
(hard gluon, hard quark, soft gluon) and the GS-HO photon parton distributions. The 
expectation from non-diffractive processes is shown as the dashed line, which is obtained 
by PYTHIA with MRSD- and GRV-HO parton distributions for the proton and the 
photon respectively. 
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Figure 50: Constraints on the J^jp — Cg plane, the momentum sum and the fraction of hard 
gluons in the pomeron. The solid line is the result of the fit of tfi® measured da/dr]^^^ 
for the diffractive event sample with the predictions of POMPYT. The two dot-dashed 
curves are obtained from the measurement of the diffractive structure function in DIS 


159|| , assuming two (lower) or three (upper) quark flavours. The shaded bands represent 


1 a contours. 
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Figure 51: r]max distribution of D* candidate events. The hatched histogram represents 
the contribution of non-diffractive processes, from PYTHIA, the solid histogram is the 
sum of the previous one with the contribution of diffractive processes from RAPGAP. 
The two components are respectively normalized to the number of events with r]max > 3 
and r]max < 2. 
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Figure 52: Resolved photoproduction via (a) colour singlet exchange and (b) colour 
non-singlet exchange. The rapidity gap event morphology is shown in (c) where black 
dots represent hnal state hadrons and the boundary illustrates the limit of the ZEUS 
acceptance. Two jets of radius R are shown, which are back to back in azimuth and 
separated by a pseudorapidity interval At]. An expectation for the behaviour of the gap 
fraction is shown in (d)(solid line). The non-singlet contribution is shown as the dotted 
line and the colour singlet contribution as the dashed line. 
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Figure 53: The gap-fraction as a function of the gap width from the 
dijet events [p.78|]. Data are compared to a £t to an exponential plus a 


ZEUS analysis of 
constant. 
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